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research  conducted  under  that  contract  from  November  1975  through  June  1976. 
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Navy  X-22A  Program  Manager  at  NADC;  in  addition,  the  authors  acknowledge  their 
appreciation  to  Mr.  C.  Mazza  at  NADC  and  Messrs.  D.  Hutchins  and  R.  Siewert  at 
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ABSTRACT 


^ A study  program  was  undertaken  to  investigate  the  feasibility  of  con- 
ducting ground-  and  in-flight  simulations  of  a jet-lift  VTOL  aircraft  performing 
terminal  area  operations  using  the  U.  S.  Navy  X-22A  V/STOL  Research  Facility. 

The  objectives  of  the  program  were: 

develop  a generic  mathematical  model  of  the  AV-8A  aircraft 
3uring  a decelerating  approach  to  hover; 

• to  develop  the  methodology  for  computing  the  variable  stability 
gains  required  to  simulate  the  AV-8A  model  with  the  X-22Ay 

f to  develop  the  procedures  for  implementing  the  simulation  in  the 
X-22A  ground  simulator© 

7 The  general  conclusion  for  this  study  is  that  it  is  feasible  to  simu- 
late in  the  X-22A  aircraft  an  AV-8A  class  VTOL  aircraft  during  terminal  area 
operations.  The  methodology  for  performing  such  a simulation  is  developed  and 
the  techniques  for  mechanization  on  the  X-22A  ground  simulator  for  further 
experimental  investigation  are  described.  Excellent  simulation  fidelity  is 
possible  during  low-speed  flight  (gj60  kts).  Areas  which  require  further  investi- 
gation include  the  discrepancies  in  response  to  throttle  inputs  at  higher  speeds 
and  the  necessity  for  matching  Lateral  acceleration  responses  at  the  pilot's 
station  to  control  inputs.  \ 
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Section  1 


INTRODUCTION 


This  report  documents  the  results  of  a study  directed  at  the  in- 
flight simulation  of  a jet-lift  VTOL  aircraft  performing  terminal  area  opera- 
tions. In  particular,  the  feasibility  of  conducting  such  a simulation  with 
the  variable  stability  X-22A  V/STOL  aircraft  was  addressed  and  the  requisite 
methodologies  were  developed.  The  study  underwent  several  changes  in  empha- 
sis and  concomitant  modification  of  the  original  objectives  during  its  perfor- 
mance period.  In  this  section,  therefore,  a brief  summary  of  the  program's 
evolution  will  be  presented  followed  by  a statement  of  the  resulting  study  ob- 
jectives; an  outline  of  this  report,  which  documents  the  accomplishment  of 
these  objectives,  concludes  the  section. 

As  originally  planned,  the  study  was  to  concentrate  primarily  on 
angular  augmentation  control  system  designs  for  a "generic"  jet-lift  aircraft; 
this  generic  model  was  intended  to  include  only  representative  values  of  major 
stability  derivatives  in  the  moment  equations,  and  hence  would  require  a mini- 
mal effort  to  formulate.  Shortly  after  the  commencement  of  the  contract, 
however,  the  effort  was  redirected  by  NAVAIR  toward  simulating  a specific  jet- 
lift  VTOL  --  the  AV-8A  Harrier  — in  order  that  the  flight  experiment  to  be 
performed  subsequent  to  this  study  would  provide  results  directly  applicable  to 
the  Harrier. 

As  a result  of  this  redirection,  considerably  more  effort  than  orig- 
inally planned  was  required  to  (1)  formulate  a suitable  Harrier  mathematical 
model  from  available  data,  and  (2)  develop  the  simulation  methodology.  The 
objectives  of  the  study  therefore  were  modified:  a major  emphasis  was  to  be 

placed  on  developing  an  accurate  model  of  the  AV-8A  so  that  one  would  be  avail- 
able at  the  onset  of  the  flight  program.  Accordingly,  for  the  first  half  of 
the  program,  attention  was  focussed  on  acquiring  and  analyzing  available  AV-8A 
data  toward  this  end. 

The  development  of  an  accurate  and  suitable  mathematical  model  was 
hampered  by  delays  in  acquiring  data  and  significant  discrepancies  among  those 
data  that  were  obtained;  further,  some  required  information  could  not  be  ex- 
tracted from  the  data  available.  In  the  same  time  period  that  the  modelling 
inadequacies  were  becoming  clear,  the  objectives  of  the  follow-on  flight  pro- 
gram were  again  altered:  rather  than  a simulation  of  the  AV-8A,  the  flight 

experiment  was  to  be  tied  into  the  development  of  the  McDonnell -Douglas  AV-8B 
aircraft,  and  the  in-flight  simulation  will  be  of  the  design  characteristics 
of  this  machine  when  they  become  available.  This  change  eliminated  the  neces- 
sity to  develop  a precise  AV-8A  model  during  this  study  program,  and  thereby 
permitted  increased  attention  to  be  paid  to  simulation  methodology  and  imple- 
mentation mechanization. 

As  can  be  seen,  the  evolution  of  the  goals  for  the  follow-on  flight 
program  led  to  concomitant  modifications  in  the  objectives  of  this  study  from 
those  given  at  its  onset.  The  final  objectives  were: 
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• Develop  a model  of  the  AV-8A  which  contains  representative 
stability  and  control  parameters  but  may  not  be  exactly  correct. 

• Develop  the  methodology  for  computing  the  variable  stability 
gains  required  to  simulate  the  AV-8A  model  with  the  X-22A. 

Assess  simulation  fidelity  to  the  extent  possible. 

• Develop  the  procedures  for  implementing  the  simulation  in  the 
X-22A  ground  simulator,  and  check  out  resulting  mechanization. 

The  remainder  of  this  summary  report  documents  the  procedures  and 
analyses  conducted  to  accomplish  these  objectives.  Section  2 discusses  the 
development  of  the  AV-8A  mathematical  model,  including  a review  of  the  requisite 
model  form  and  the  form  in  which  the  available  data  appeared,  a summary  of 
necessary  assumptions  and  procedures,  and  a documentation  of  the  resulting  model 
characteristics.  A general  discussion  of  means  for  investigating  and  assessing 
simulation  fidelity  is  presented  in  Section  3.  In  Section  4,  detailed 
descriptions  of  the  development  and  application  of  various  methods  for  deter- 
mining the  variable  stability  system  gains  for  this  problem  are  given,  both  for 
the  lateral-directional  simulation  (4.1)  and  for  the  longitudinal  computations 
(4.2).  The  modifications  made  to  the  ground  simulator  and  the  implementation 
of  the  AV-8A  model  simulation  are  described  in  Section  5,  including  the  check- 
out results.  Finally,  Section  6 summarizes  the  conclusions  of  this  study  and 
the  implications  of  the  results  for  the  follow-on  AV-8B  in-flight  simulation. 
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Section  2 


DEVELOPMENT  OF  AV-8A  MATHEMATICAL  MODEL 


2.1  REQUIRED  MATH  MODEL  AND  EXISTING  DATA  FORM 

The  variable  stability  system  of  the  X-22A  is  of  the  response  feed- 
back type.  As  implemented  in  the  X-22A,  the  characteristics  of  such  a system 
imply  certain  constraints  on  the  model  of  the  aircraft  to  be  simulated.  In 
particular,  the  system  commands  designated  control  motions  linearly  proportional 
to  measured  signals  less  "reference"  values  (which  may  be  constant  or  a function 
of  flight  condition)  to  alter  the  basic  stability/control  derivatives  of  the 
X-22A  directly;  the  proportionality  between  a measurement  and  a control  motion 
may  be  varied  as  a function  of  flight  condition  but  it  is  linear  (e.g.,  the 
control  cannot  move  as  the  square  of  the  measurement).  As  a result,  the  model 
of  the  AV-8A  must  be  developed  in  terms  of  reference  conditions  which  explicitly 
depend  on  a defined  evaluation  task,  and  "linearized"  (though  non-constant) 
stability  and  control  derivatives  around  these  references. 

As  an  example,  consider  the  longitudinal  equations  of  motion  which, 
with  some  simplifications,  may  be  written: 

tc  ■*  q Sen  & *•  Our  - X ( u. , US , q , ST  , S g , Gj  ) 

ulr  - q COS  Q - G CL  - 2 ( OL,  -ocr,  q , St  , Se  , £j) 

6 - rt(u.,  CO,  q , Sr  , <$e  . dj) 

y s?  0 - at 

at  - ta.n  ' ( -ccr / ol  ) 

By  requiring  that  0 = G=  <p  = y = Ot*  Obe  included  as  reference  conditions, 
this  set  of  equations  may  be  reduced  to  three  equations  in  seven  unknowns: 

U > y , U , or , , St  , Gj  . By  specifying  four  of  the  unknowns  through 

a selection  of  the  evaluation  task,  the  -emaining  three  may  be  uniquely  deter- 
mined. For  example,  the  reference  conditions  at  various  points  along  the  AV-8A 
reference  trajectory  may  be  computed  by  specifying  y , at  , u and  6j  and  cal- 
culating the  required  values  of  ii , Sr  , and  c$g  from  the  complete  "global" 
representation  of  the  aerodynamic/propulsive  forces  and  moments  acting  on  the 
aircraft.  If  the  deceleration  schedule  is  defined  for  the  task,  then  u.  , y , 
a and  Oj  (or  aat)  may  be  specified  and  the  reference  values  of  ST  , cft  and 
*tr  (or  &j  ) may  be  calculated.  These  reference  conditions,  either  specified 
or  calculated,  are  denoted  by  the  subscript  R for  convenience. 
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The  coefficients  for  a Taylor  series  representation  of  the  aerodynamic/ 
propulsive  forces  and  moments  may  now  be  obtained  in  a manner  similar  to  the 
technique  described  in  Reference  1.  The  general  approach  used  to  calculate  the 
longitudinal  dimensional  derivatives  is  to: 

• calculate  the  values  of  X , Z , and  M which  correspond  to  the 

reference  conditions  [ a .R  , uJg  , (f  R , Srg  , 6eg  > at 

selected  points  on  the  reference  trajectory 

• perturb  one  of  the  variables  in  suitable  increments  from  the 
reference  value  and  determine  the  corresponding  values  of 

X , Z and  A 1 . 

• calculate  and  plot  the  change  in  X,Z  and  versus  the  change 
in  the  variable  being  perturbed 

• fit  an  n^  order  polynomial  to  the  data  on  a least-square 
error  basis  (e.g.  n = 3) 

• differentiate  the  polynomial  successively  to  determine  first, 
second,  and  third-order  derivatives. 

The  cross  derivatives  such  as  *t*jrmay  be  obtained  by  repeating  the  above 
procedure  for  a series  of  throttle  positions  at  positive  and  negative  increments 
from  the  reference  value,  cross-plotting  the  resultant  values  of  A^vs.  ST  , 
and  calculating  the  slope  of  the  curve. 

It  is  clear,  then,  that  the  development  of  an  AV-8A  model  in  the  appro- 
priate form  is  an  involved  process  unless  data  in  a similar  form  are  available 
The  data  reports  investigated  for  this  study  are  listed  below  (References  2 
through  21,  respectively): 

MCAIR  A1410:  AV-8A  Aerodynamics  (Wind  Tunnel  Data  Summary) 

MCAIR  A1411:  AV-8A  Estimated  Flying  Qualities 
Rockwell  NR72H-268:  AV-8A  Hover  Simulation 
MCAIR  A3278:  AV-8A  Landing  Approach  Simulation 
MCAIR  A3618:  AV-8A  Landing  Approach  Simulation 
MCAIR  A3S80:  AV-8A  Stability  Augmentation  Design 
NRC  10257:  XV-6A  Landing  Approach  Simulation 
NASA  TND-6791:  XV-6A  Landing  Approach 
AFFTC  TR  No.  69-26:  USAF  Harrier  Evaluation 
NATC  FT-98R-66 : XV- 6 A Carrier  Suitability 
AFFTC  TR  No.  68-iO:  XV- 6 A Handling  Qualities 
NATC  FT-52R-71:  Harrier  Flying  Qualities 
NATC  FT-59R-71:  AV-8A  Carrier  Suitability 

NASA  TND-7296:  XV-6A  Aerodynamic  Parameters  from  Flight  Data 

NASA  TND-6826 : XV-6A  Wind  Tunnel  Study 

MCAIR  A3091:  AV-8A  Control  System  Analysis 

NATC  FT-72R-72:  AV-8A  Shipboard  Suitability 

NATC  FT-05R-69 : Navy  MPE  of  Harrier 

NATC  FT-27R-71 : AV-8A  Flying  Qualities 

NAVAIR  01-AV8A-1 : AV-8A  Flight  Manual 
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References  4-8  contain  scattered  documentation  on  stability  and  control 
derivatives  in  the  desired  dimensional  form  (e.g.,  etc.)  as  a function  of 
velocity;  these  data  do  not  include  trim  information  and  are  generally  incon- 
sistent with  each  other.  The  flight  test  reports  (References  9 - 14,  18  - 20) 
are  useful  for  qualitative  understanding  of  trim  requirements  and  transition 
capabilities,  although  again  inconsistencies  appear  in  the  limited  quantitative 
information  that  exists.  Reference  2 is  a compilation  of  AV-8A  wind  tunnel  data 
in  non-dimensional  form,  generally  for  speeds  above  50  kts ; it  is  the  most 
complete  of  the  data  sets  available,  although  considerable  "massaging"  is  re- 
quired to  obtain  the  desired  form,  and  it  gives  a qualitative  understanding  of 
the  global  characteristics  of  the  AV-8A.  Reference  3 presents  characteristic 
roots  computed  on  the  basis  of  the  Reference  2 data;  although  incomplete,  some 
comparisons  of  the  stability  characteristics  of  a developed  model  can  be  made. 

A general  summary  of  the  quality  and  deficiencies  of  the  available 
data  is  therefore: 

• Stability  derivatives.  Fairly  complete  information  from 
Reference  2 above  50  kts.  Comparison  of  selected  derivatives 
possible  with  References  7,  8;  comparison  of  selected  modal 
characteristics  possible  with  Reference  3. 

• Control  derivatives.  Information  about  the  three  moment  con- 
troller effectivenesses  available  in  Reference  2.  Essentially 
no  quantitative  information  on  throttle  sensitivities  for  the 
longitudinal  force  equations  or  on  nozzle  angle  effectiveness; 
no  information  on  moment  sensitivities  for  thrust  and  nozzle. 

• Reference  values.  No  direct  information.  References  5 and  6 
discuss  landing  approach  evaluation  tasks  from  which  some 
reference  information  may  be  inferred;  limited  quantitative 
data  from  flight  test  are  in  References  10,  12,  18  - 20. 

• Trims.  No  direct  information  for  throttle.  Longitudinal 
stick  trim  estimates  may  be  inferred  from  Reference  3. 

2.2  AV-8A  MATH  MODEL  GENERATION 

As  the  preceding  discussion  indicates,  significant  gaps  exist  in  the 
data  base  available  for  formulating  an  AV-8A  mathematical  model  in  an  appropriate 
form.  In  view  of  the  redirection  to  the  AV-8B  aircraft  for  the  follow-on  flight 
program,  it  was  decided  to  formulate  the  AV-8A  model  for  this  study  by  extra- 
polating from  these  data  where  possible  and  estimating  undocumented  character- 
istics by  "engineering  judgment"  (e.g.,  geometries).  A summary  of  the  resulting 
assumptions  and  procedures  is  given  in  this  subsection. 

The  initial  assumptions  involve  the  definition  of  an  evaluation  task, 
since  the  reference  conditions  depend  upon  it.  Following  the  general  guidelines 
provided  by  previous  MCAIR  simulations  (References  5,  6),  the  following  task 
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characteristics  were  selected: 


• constant  glide  slope:  ~f  = -5  deg 

• constant  angle  of  attack:  (X  - 8 deg 

• two  possible  initial  velocities:  VQ  = 80  kt,  120  kt 

• one-step  deceleration:  nozzle  angle  ^ 81  deg  "instantaneously" 
at  appropriate  range 

The  choice  of  8 degrees  constant  angle  of  attack  rather  than  pitch  attitude  was 
selected  by  MCAIR  because  the  AV-8A  runs  out  of  longitudinal  control  authority 
at  approximately  15  degrees  for  low  speeds;  the  two  approach  speeds  are  those 
investigated  by  MCAIR  in  Reference  5. 

Longitudinal  and  lateral-directional  stability  derivatives  were  then 
calculated  from  the  non-dimensional  wind  tunnel  data  of  Reference  2 for  six 
trim  (unaccelerated)  flight  conditions,  all  at  f = -5  deg,  ct  = +8  deg:  VQ  - 0, 
30,  60,  80,  100,  120  kt.  The  resulting  dimensional  stability  derivatives  are 
given  in  Appendix  I;  also  shown  are  available  data  from  References  4,  6,  7, 
and  8 when  applicable.  It  is  emphasized  that  the  V0  - 0,  30  kt  data  are  extra- 
polations: the  wind  tunnel  data  in  Reference  2 are  generally  valid  only  for 

/ e 60  kt.  In  general,  consistent  trends  are  exhibited,  and  the  agreement 
among  data  sets  is  fair  (with  the  exception  of  some  Reference  4 data  incorrect 
signs).  Unfortunately,  the  major  uncertainties  are  for  the  derivatives  M u , 

^v- » 311(1  *v  • which  have  a major  effect  on  the  aircraft  dynamics;  the 
selectea  values  of  *40  and  are  a reasonable  average  of  the  Reference  2 and 
Reference  8 values,  while  the  Lv  and  A/^,  values  of  Reference  7 were  weighted 
most  highly  in  the  selection  of  these  derivatives. 

Table  2-1  summarizes  the  stability  derivatives  of  the  resulting  AV-8A 
model  used  in  this  study.  Figure  2-1  compares  the  characteristic  roots  of  this 
model  at  the  six  trim  velocities  with  available  documentation  from  Reference  3; 
considering  the  assumptions  and  extrapolations  made  in  the  model  development, 
the  correlation  is  outstanding  for  both  the  longitudinal  and  the  lateral- 
directional  degrees  of  freedom.  It  appears  that  this  part  of  the  developed 
AV-8A  model  is  representative  of  the  MCAIR  estimate  of  the  actual  aircraft. 

The  moment  control  derivatives  were  also  estimated  from  the  data  pre- 
sented in  Reference  2.  Both  aerodynamic  and  reaction  control  data  are  required: 
the  reaction  control  derivatives  were  selected  for  maximum  thrust  and  singular 
demand  (no  effect  of  simultaneous  control  usage  on  available  bleed  air  is 
included).  Appendix  I presents  the  data;  Table  2-2  summarizes  the  estimated 
control  effectivenesses.  The  values  are  somewhat  smaller  than  those  given  in 
Reference  7,  and  hence  may  be  considered  to  reflect  some  effects  of  simultaneous 
control  usage. 

Quantitative  throttle  control  effectiveness  data  are  very  limited  in 
the  available  references.  The  value  of  Fg  for  the  nozzles  at  90  degrees  was 
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TABLE  2-1 

AV-8A  MODEL  STABILITY  DERIVATIVES 


V 

Derivative 

0 

30 

60 

80 

100 

120 

- .026 

— 

.035 

- .043 

— 
- .049 

- .055 

- .061 

0.0 

- 

.020 

- .027 

- .029 

- .030 

- .030 

0.0 

- 

.011 

- .036 

- .060 

- .083 

- .108 

- .020 

- 

.105 

- .190 

- .250 

- .310 

- .370 

*u 

.001 

.0017 

.0022 

.0025 

.0026 

.0026 

.0009 

.0042 

.0054 

.0059 

.0060 

.0058 

% 

0.0 

- 

.13 

- .26 

- .34 

- .43 

- .52 

V 

- .025 

" 

.042 

- .059 

- .070 

- .081 

- .092 

C 

- .005 

- 

.012 

- .04 

- .053 

- .059 

- .058 

4 

- .04 

- 

.26 

- .48 

- .64 

- .80 

- .94 

c 

0.0 

.110 

.215 

.285 

.360 

.430 

< 

- .0045 

- 

.0021 

0.00 

.0013 

.0019 

.0020 

< 

0.0 

- 

.075 

- .145 

- .195 

- .240 

- .290 

< 

- .015 

- 

.060 

- .100 

- .130 

- .160 

- .190 
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TABLE  2-2 

AV-8A  MODEL  MOMENT  CONTROL  DERIVATIVES 


V 

0 

Derivative* 

30 

60 

80 

100 

120 

IU 

*o 

X 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

?<5e 

.41 

.53 

.67 

.77 

.88 

1.00 

f«8e 

.15 

.20 

.25 

.28 

.31 

.34 

.45 

.45 

.45 

.45 

.45 

.45 

.014 

.014 

.014 

.014 

.014 

.014 

- .61 

- .69 

- .78 

- .84 

- .90 

- .96 

LS* 

- .080 

- .085 

- .095 

- .100 

- .105 

- .110 

.185 

.215 

.240 

.260 

.280 

.295 

•per  inch  control  deflection 
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selected  to  be  0.1  g/in.  =3.2  ft/sec2/in.  This  value  is  higher  than  that  used 
in  Reference  4 (2.58),  Reference  8 (2.7),  and  an  estimate  based  on  data  from 
Reference  9 (2.6),  but  is  lower  than  estimates  made  from  information  given  in 
Reference  14  (4.9)  and  Reference  21  (5.6);  no  direct  information  on  the  thrust 
derivatives  is  given  in  the  MCAIR  simulation  reports  (References  5,  6,  7).  With 
this  value  selected,  the  values  of  and  at  other  nozzle  angles  were 
estimated  by  simple  geometric  arguments: 

- ~ 3.2  [cos  (90  - 0j)] 

XST  ~ s-  2 [sen  o°  ~ 0j  )] 

/ Q 

where  Gj  = + 2 (engine  angle  in  fuselage) 

Essentially  no  information  on  pitching  moment  caused  by  thrust  was  available 
in  the  literature.  Although  a value  of  = -0.03  rad/sec2/in.  could  be  com- 
puted from  Reference  2 for  one  flight  condition,  the  sign  of  this  estimate  was 
inconsistent  with  qualitative  information  from  flight  test.  It  was  therefore 
assumed  arbitrarily  that  M$r  was  positive,  increased  with  speed,  and  decreased 
with  increasing  nozzle  angle.  The  resulting  values  of  X&  , it  , and  Hj  are 
listed  in  Table  2-3.  77  7 

The  remaining  formulation  required  is  the  trim  information  --  the 
reference  values  for  the  controls.  No  antisymmetric  characteristics  of  the  air- 
craft are  included  in  the  model,  and  so  the  trim  lateral  stick  and  rudder  pedal 
control  positions  are  zero.  Longitudinal  stick  trim  data  were  available  in 
Reference  3,  and  are  reproduced  in  Figure  2-2.  Also  shown  is  an  approximate 
variation  labelled  - 0 calculated  on  a piece-wise  linear  basis  using  the  form- 
ulated stability  and  control  derivatives  and  reference  conditions  from  the  equation: 

AS  - ~ Mu  Aa  - a ur 

£ 

where  = “trim,  ~ ’ etC> 

As  can  be  seen,  the  correct  trim  trend  is  reproduced  although  somewhat  larger 
changes  are  predicted;  the  differences  between  the  curves  can  be  used  as  a 
qualitative  check  on  the  values  of  assumed  earlier  once  the  trim  throttle 
positions  are  computed. 

The  trim  throttle  positions  and  nozzle  angles  were  estimated  on  the 
basis  of  geometric  arguments  since  again  no  data  were  available  in  the  refer- 
ences. On  a piece-wise  linear  basis,  the  following  approximate  expressions 
were  used  to  compute  the  required  X and  Z forces  supplied  by  the  propulsion 
system: 
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where 


AX  = 


AZ  = 


Zu  Au  + 

A ( 


Xtxr  bu/  - b(iio  -b  g sin  G0) 

2^  A*"  + Zses  &6e  ~ b(to0  -g  cos  0O) 


Approximate  values  for  the  kinematic  terms  were  obtained  by  using  the  reference 
condition  cc0  = 8 degrees,  Y = -5  degrees  and  assuming  an  average  deceleration 
of  approximately  0.05  g on  the  glide  slope.  From  these  computations,  a plot 
of  X vs.£*  propulsive  force  is  shown  in  Figure  2-3.  Using  6*y  = 81  degrees  as 
whe  hover  nozzle  angle,  approximate  trim  nozzle  angles  for  the  velocities  of 
interest  can  be  found.  With  these  values,  trim  throttle  positions  may  be  computed 
using  the  relationships  for  XgT  or  lg  given  earlier;  the  approximate  trim  throttle 
positions  are  plotted  in  Figure  2-4. 


As  was  noted  earlier,  it  is  possible  to  check  qualitatively  the 
assumed  values  of  Mgr  by  using  the  developed  trim  throttle  and  longitudinal 
stick  positions  along  with  the  estimate  of  longitudinal  control  sensitivity: 


Me  A &e 

~ a£  er«/M 

^ r TRIM 

Table  2-4  presents  the  comparison  of  these  values  with  the  assumed  values  of 
Mgr  , and  indicates  qualitative  agreement  in  magnitude  and  sign.  Hence,  the 
assumed  values  are  at  least  reasonable,  and  were  therefore  retained  in  the  model. 


The  formulation  of  a generic  AV-8A  model  is  thus  complete.  Table 
2-5  summarizes  the  transfer  functions  of  the  model  at  three  flight  conditions 
for  reference.  The  simulation  of  this  model  with  the  X-22A  is  discussed  in 
Sections  4 and  5 of  this  report. 
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ft/sec 
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TABLE  2-4 

COMPARISON  OF  ASSUMED  M,$r  ( 03  = 81°) 

WITH  APPROXIMATE  VALUE  DERIVED  FROM  PITCH  STICK  TRIM 


V/ 

Assumed 

Derived 

0 

0.0 

0.0 

30 

.005 

.006 

60 

.010 

.014 

80 

.015 

.025 

100 

.020 

.020 

120 

.020 

0.0 

TABLE  2-5 

TRANSFER  FUNCTIONS  OF  AV-8A 
MODEL  IN  FORM  K (f/r)  [<". 


■■ 

V - 0 kt*  . 

V = 60  kt* 

V = 120  kt** 

oho 

(.33) (.02)  [-.48;.3l3 

(1.00)(.073).Q-.91;.32] 

(1.53)  (.084) (- .46) (- . 187) 

-4.8(.022) 

-3. 6(2. 5) (.19) 

-9.69(.31) (1.42) 

"Z 

. 41  (. 33)  [-.48;.3l] 

.67(37.9)  Q.14;.ll] 

1.0(68.5)  Q 26;. 13^ 

. 15 ( . 022) (.026) (0) 

. 25 ( . 21) (.036) (0) 

. 34 (.  40) ( . 051) (0) 

.38(.63)  [-  .49;  .62^ 

.38(1.45)  Q-.61;.95] 

1.97(1.51) (-.87) (-.075) 

< 

-3. 18(.33)  [-.48;.3l] 

-3 . 18 ( . 40)  [-.S0;.4l] 

-2. 52 (.33) (-1.55) (-.10) 

N*r 

-.0025 (.027) (0) 

. O' ( . 027) (-1.43) (0) 

.025  [.19;. 13]  (0) 

A L/D 

(.015)(.58)[-.S0;.52] 

( .068) (1 . 26)  [-.28;1.2j 

(.056) (1.73)  [-.15;1.8] 

nT 

. 4 ( .015) (34 . 3) 

5 .0 ( .075) (4 . 15) 

10.0(.074) (3.99) 

< 

.45(-.17)  (.21) 

.45  [.23;.  36]] 

.45  [.17;. 883 

.014[.53;1.64j(-1.66) 

.014(.49)(-.66)(-4.0) 

,014(.53)(-.94)(-7.9) 

- .6(.015)  £.67;2.07j 

-.78(.08)(.69)(32.8) 

- . 96 (- .052) (1.06) (64 . 8) 

<r 

-.0o(.37) (-.37) 

- . 10(2.9) (-3.6) 

- . 11 (4 .9) (-6.3) 

"*rp 

.19  [-.46;.6j  (.63) 

.24  [- . 24 ; 1.07^  (1.11) 

.3  [.01 ; 1 . 33  (1.12) 

* 0J  = 81° 

**  * 50° 
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Section  3 


SIMULATION  FIDELITY 


3.1  INTRODUCTION 

The  purpose  of  this  section  is  to  present  a discussion  of,  and  a 
rationale  for,  the  guidelines  used  to  establish  the  preliminary  criteria  for 
a valid  piloted  simulation  of  the  AV-8A  aircraft  during  a decelerating  instru- 
ment landing  approach.  "It  is  trite  to  say  that  to  build  an  exact  replica  of 
the  aircraft  to  be  simulated,  and  operate  it  in  a representative  environment, 
is  a perfect  simulation,  and  yet  no  simulation  at  all."  (Reference  22).  This 
statement  contains  the  primary  issue  of  the  subject  of  simulation  fidelity: 
since  exact  replication  is  not  possible,  a determination  must  be  made  of  which 
characteristics  one  should  attempt  to  duplicate,  and,  further,  to  what  accuracy 
the  duplication  must  be  performed.  The  extent  to  which  this  determination  can 
be  objective  rather  than  subjective  is  dictated  by  the  exigencies  of  the  par- 
ticular simulation  problem:  imperfect  knowledge  about  the  effects  of  imperfect 

duplication  of  all  simulation  characteristics  renders  a generalized  set  of  ob- 
jective measures  impossible.  This  inchoate  situation  has  a direct  parallel  - 
and  raison  d'etre  - in  the  continuing  research  efforts  to  devise  appropriate 
parameters  for  flying  qualities  specifications,  and  fundamentally  mandates  a 
determination  process  that  is  incompletely  objective  and  at  least  partly 
judgmental.  Factors  which  have  a major  influence  on  the  objectivity  possible 
in  assessing  simulation  validity  include: 

• Objectives  of  the  simulation 

training/procedures  vs.  research/flying  qualities 
missions,  tasks 

• Quality  of  AV-8A  mathematical  model 

validity  of  model  used 
accuracy  of  individual  parameters 

• Simulator  capabilities  and  limitations 

ability  to  reproduce  model  stability/control  characteristics 
motion/visual  cues 
ergonomic  considerations 
environmental  factors. 

The  implications  of  the  applicable  topics  subsumed  by  each  of  the 
above  three  categories  for  the  AV-8  simulation  problem  are  briefly  discussed 
in  the  following  subsections;  following  this  discussion,  the  resultant  general 
guidelines  from  which  were  derived  the  specific  techniques  for  implementing  the 
AV-8A  model  on  the  X-22A  simulator/aircraft  are  summarized. 
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3.2 


OBJECTIVES  OF  THE  SIMULATION 


The  general  goal  of  the  in-flight  simulation  experiment  for  which  this 
program  is  a prerequisite  is  to  provide  information  relative  to  the  interaction 
of  angular  augmentation  control  systems,  displayed  information,  and  task 
variables  during  the  landing  approach  of  a jet-lift  VTOL  aircraft. 

General  guidance  for  the  accuracy  of  VTOL  aircraft  flying  qualities 
simulations  such  as  the  subject  AV-8A  simulation  is  implicitly  contained  in 
the  military  flying  qualities  specification,  MIL-F-83300  (Reference  23  ),  and 
suggested  revisions  (Reference  24).  These  documents  identify  key  flying 
qualities  parameters  for  both  forward  and  hover  flight  and  delineate  boundaries 
for  these  parameters  as  a function  of  aircraft  Class  and  Flight  Phase  Category 
based  upon  pilot  rating  data.  This  concept  of  flying  qualities  Levels  may  also 
be  applied  to  the  simulation  fidelity  problem;  the  important  flying  qualities 
parameters  of  the  AV-8A  model  and  the  X-22A  simulation  of  that  model  must  be 
within  the  same  flying  qualities  Level  boundaries  in  order  that  a valid  simula- 
tion can  reasonably  be  anticipated. 

A similar  technique  has  been  applied  to  the  task  of  establishing 
dynamic  response  criteria  for  simulators  used  in  training  applications  in  Refer- 
ence 25  for  conventional  aircraft  and  in  Reference  26  for  VTOL  aircraft. 
Allowable  tolerances  in  significant  flying  qualities  parameters  are  established 
based  upon  pilot  rating  sensitivity  analyses;  the  range  of  acceptable  values 
for  each  parameter  corresponds  to  that  interval  of  values  about  the  nominal  for 
which  pilot  rating  varies  no  more  than  one  rating  unit  (Reference  25  ) , or  one- 
half  a rating  unit  (Reference  26  ) , about  the  nominal  rating.  Although  flying 
qualities  research  programs  in  general  require  a higher  degree  of  simulation 
fidelity  than  training  simulations,  the  results  of  Reference  26  as  well  as 
MIL-F-83300  provide  useful  guidance  for  the  AV-8A  simulation  problem  in  cases 
where  a particular  stability/control  parameter  of  the  AV-8A  cannot  be  exactly 
reproduced  by  the  X-22A  because  of  system  limitations. 

The  specific  objective  of  the  flight  program  is  to  define  the  dynamic 
characteristics  required  of  the  various  stability/control  augmentation  systems 
as  a function  of  major  task  elements  such  as  breakout  range/altitude  for  rep- 
resentative levels  of  displayed  information  (probably  both  head-up  and  head- 
down).  Although  the  evaluation  task  details  are  currently  undefined,  it  is 
reasonable  to  assume  an  approach  profile  of  the  type  used  in  MCAIR  ground  simu- 
lations (Reference  6 ) , which  consist  of  a level  deceleration  to  partially 
jet-borne  flight,  constant  speed  instrument  glide  slope  tracking,  and  a final 
deceleration  just  prior  to  or  during  breakout.  The  purpose  of  the  augmentation 
systems  to  be  considered  is  to  improve  the  "short-term"  attitude  response 
characteristics,  and  the  evaluation  procedure  will  emphasize  tracking  tasks 
which  highlight  the  short-term  behavior.  Hence,  short-term  aspects  of  the 
response  to  control  inputs  (e.g.  the  first  five  seconds)  should  be  matched  as 
accurately  as  possible  with  less  weight  attached  to  longer-term  matching. 
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Stability  problems  wnich  are  inherent  in  the  jet-lift  VTOL  concept 
must  be  accurately  simulated.  The  we  11 -documented  lateral-directional  flying 
qualities  deficiencies  of  the  Harrier  between  90  and  30  knots  (Reference  21  , 

for  example)  caused  by  large  stable  values  of  L @ and  small  values  of 
(unstable  below  ~60  knots)  must  be  very  nearly  duplicated  in  the  X-22A  simu- 
lation. Longitudinally,  the  AV-8A  exhibits  an  angle-of-attack  instability 
throughout  the  transition  speed  regime  which  must  also  be  reproduced  in  the 
simulation. 

The  foregoing  considerations  have  all  been  pertinent  to  a fixed 
operating  point  simulation  of  the  Harrier.  However,  a primary  part  of  the 
evaluation  task  involves  a decelerating  transition  during  the  landing  approach; 
hence  the  important  characteristics  of  the  AV-8A  during  a conversion  from  par- 
tially jet-borne  flight  to  jet-borne  flight  must  also  be  simulated.  Experi- 
ence in  previous  X-22A  programs  (References  1 and  27  ) has  shown  that  transi- 
tions at  relatively  low  deceleration  levels  may  be  very  adequately  treated  as 
a series  of  fixed  operating  points  defined  by  a single  variable,  e.g.  X-22A 
duct  angle,  and  that  no  explicit  deceleration  effects  need  to  be  included  in 
the  equations  of  motion.  Two  characteristics  of  the  Harrier  in  transition 
must,  however,  be  included  in  the  simulation: 

• Aircraft  response  to  nozzle  rotation.  The  primary  AV-8A  deceler- 
ation control  is  through  thrust  vector  rotation.  The  celerity 
with  which  the  AV-8A  nozzle  angle  may  be  varied  requires  a con- 
sideration of  simulating  not  only  the  long-term  deceleration  char- 
acteristics of  the  vehicle,  dictated  primarily  by  the  AV-8A  drag 
damping,  but  also  the  short-term  deceleration  caused  by  the  almost 
instantaneous  rotation  of  the  thrust  vector. 

• Trim  changes  in  transition.  It  is  desirable  to  reproduce  the  gen- 
eral characteristics  of  pitch  stick  and  throttle  trim  changes  with 
speed  and  nozzle  setting  in  transition  in  order  that  the  control 
task  be  similar  to  the  AV-8A.  Perfect  replication  is  less  impor- 
tant than  for  the  constant  speed  characteristics,  however,  as  long 
as  the  general  sense  of  the  motion  required  is  correct,  because 
the  transitions  are  non -constant  maneuvers. 

3.3  QUALITY  OF  THE  AV-8A  MATHEMATICAL  MODEL 

The  concept  of  using  a simulation  technique  based  upon  linearized, small 
perturbation  equations  of  motion  for  a VTOL  aircraft  in  transition  provides  an 
adequate  ground  simulation  of  the  X-22A  (Reference  27).  However,  a valid  simu- 
lation of  the  Harrier  in  transition  may  need  to  include  some  important  higher- 
order  effects;  for  example: 

• The  effect  of  thrust  level  on  control  authority/sensitivity.  As 
indicated  in  Appendix  I,  RCS  control  power  is  dependent  upon  thrust 
level.  Maximum  thrust  has  been  assumed  in  the  calculation  of  the 
model  control  derivatives;  the  effect  on  control  power  of  any  large 
deviation  from  maximum  thrust  may  have  to  be  accounted  for  in  the 
simulation. 


• The  effect  of  multi-control  usage  on  individual  control  authori- 
ties/sensitivities. Appendix  I also  demonstrates  the  reduction 

in  individual  control  power/sensitivity  which  results  from  a simul- 
taneous usage  of  one  or  both  of  the  remaining  RCS  controls.  The 
inclusion  of  this  effect  may  be  required  in  order  to  ensure  simu- 
lation fidelity  in  the  event  of  large  deviations  from  "trim"  con- 
trol positions. 

• The  effect  of  large  perturbations  of  the  state  variables.  The 
sensitivity  of  any  of  the  stability/control  parameters  of  the 
AV-8A  presented  in  Section  2 to  perturbations  from  the  established 
reference  conditions  is  currently  unknown.  Conversations  with 
MCAIR  personnel  have  indicated  that,  for  example,  the  value  of 

is  heavily  dependent  upon  angle-of-attack  and  power  setting. 
These  effects  cannot  currently  be  quantified  and  evaluated;  hence 
they  are  not  included  in  the  simulation.  However,  future  simula- 
tions of  the  AV-8B  will  require  an  investigation  of  these  effects. 

In  addition  to  the  exclusion  of  the  above  higher-order  effects,  a fur- 
ther potential  limitation  to  the  fidelity  of  the  simulation  is  the  accuracy  of 
the  AV-8A  model  stability/control  parameter  values.  As  discussed  in  Section  2, 
these  values  are  the  result  of  a fairly  extensive  data  search,  analyses,  and  the 
application  of  engineering  judgment;  the  object  is  to  provide  the  best  possible 
simulation  of  the  resulting  AV-8A  model  characteristics  rather  than  to  be  con- 
cerned about  the  detailed  validity  of  the  model. 

3.4  SIMULATOR  CAPABILITIES  AND  LIMITATIONS 

The  X-22A  aircraft,  hence  simulator,  does  not  have  independent  con- 
trol over  all  six  degrees  of  freedom.  This  fundamental  limitation  has  impli- 
cations for  both  the  longitudinal  and  lateral-directional  simulation  of  the 
AV-8A.  Longitudinally,  the  X-  and  Z-force  characteristics  of  the  AV-8A  cannot 
be  independently  matched  by  the  X-22A  since  collective  blade  pitch  is  currently 
the  only  longitudinal  force  control  available,  and  its  effect  is  dictated  pri- 
marily by  thrust  vector  (duct)  angle.  The  ability  to  incorporate  the  collec- 
tive use  of  the  elevons  as  a third  longitudinal  controller  will  serve  to  amel- 
iorate this  potential  problem  area;  however,  no  means  exists  for  a similar 
alleviation  in  the  lateral-directional  degrees  of  freedom.  The  higher  value 
of  side-force-due-to-sideslip  of  the  X-22A,  which  is  not  variable,  means  that 
matching  both  sideslip  and  lateral  acceleration  time  histories  perfectly  is  not 
possible  if  bank  angle  responses  are  matched.  The  effect  of  a lateral  accel- 
eration mismatch  will  not  be  revealed  in  a fixed-base  ground  simulator  experi- 
ment; in  flight,  lateral  acceleration  may  be  an  important  proprioceptive  and 
vestibular  sideslip  cue,  and  continuous  control  of  sideslip  is  important  in  the 
AV-8A  at  low  speeds  because  of  the  high  dihedral  effect  and  reduced  roll  control 
power.  As  a result,  a compromise  between  matching  lateral  acceleration  and 
sideslip  response  may  be  required  for  a valid  in-flight  simulation. 

The  effect  of  the  discrepancy  between  the  rotation  rates  of  X-22A 
duct  angle  and  AV-8A  nozzle  angle  has  been  discussed  as  it  applies  to  the  sim- 
ulation of  the  initial  deceleration  due  to  an  instantaneous  nozzle  rotation. 
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In  addition,  the  existence  of  a variable  mismatch  between  AV-8A  nozzle  angle 
(#y  ) and  X-22A  duct  angle  (A)  through  a transition  may  require  that  the 
simulator  throttle  gearings  be  made  functions  of  the  difference  between#,  and 
> rather  than  simply  functions  of  A . 

In  addition  to  the  well -documented  capabilities  of  the  X-22A  VSS  and 
the  additional  flexibility  afforded  by  the  airborne  analog  computer,  the  simula- 
tor and  aircraft  are  both  equipped  with  a variable  feel  system  which  adds  an 
important  dimension  to  the  potential  fidelity  of  the  AV-8A  simulation.  A faith- 
ful reproduction  of  the  AV-8A  control  system  force-feel  characteristics  in  pitch, 
roll,  and  yaw  is  possible,  given  sufficient  information  about  them.  Moreover, 
the  capability  to  vary  the  control  trim  rates,  especially  in  pitch,  yields  the 
ability  to  duplicate  the  Harrier  trim  control  characteristics  during  the  tran- 
sition. 

3.5  SIMULATION  GUIDELINES 

Based  upon  the  above  considerations  of  which  AV-8A  characteristics 
should  be  reproduced  for  a valid  simulation  and  which  of  these  desiderata  can 
reasonably  be  accomplished  within  the  current  limitations  of  the  X-22A  research 
facility,  the  following  guidelines  for  simulation  fidelity  are  proposed: 

• Fixed  operating  point  simulation  (V  = 0,  30,  60,  80,  100,  120  kts) 

1.  Match  control  system  force-displacement  relationships. 

Match  control  authorities. 

2.  Match  short-term  ( ~5  seconds)  pitch,  roll,  and  yaw  atti- 
tude reponses  of  the  simulator  and  the  AV-8A  mathematical 
model  to  specified  control  inputs. 

3.  Match  major  modal  characteristics  (short  term  frequency, 
damping;  Dutch  roll  frequency,  damping;  roll  mode  time  con- 
stant; roll-due- to-sideslip) . 

4.  Attempt  to  reproduce  specific  derivatives  judged  to  be  im- 
portant, e.g.  ^ , /V^  . 

5.  Attempt  to  match  short-term  vertical  and  longitudinal  velo- 
city responses  to  throttle  inputs,  i.e.,  time  constants  and 
relative  magnitudes. 

6.  Attempt  to  match  short-term  lateral  velocity  response  to 
control  inputs.  Further  study  of  the  lateral  acceleration/ 
sideslip  mismatch  is  warranted. 

• Transition  Simulation 

1.  Treat  transition  as  series  of  fixed  operating  points  where 
dynamic  response  characteristics  are  functions  of  duct 
angle,  hence  trimmed  velocity.  Consider  treating  throttle 
gearings  as  functions  of  (fly -A). 
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2.  Match  trends  of  trim  changes  (pitch,  power)  in  transition. 

3.  Attempt  to  match  both  short-term  and  long-term  deceleration 
characteristics  in  response  to  nozzle  rotation. 

The  effects  of  these  guidelines  on  generating  a simulation  methodol- 
ogy and  on  determining  the  techniques  for  the  actual  mechanization  of  the  AV-8A 
simulation  are  described  in  the  following  two  sections. 
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Section  4 


SIMULATION  METHODOLOGY 


4.1  LATERAL-DIRECTIONAL 

As  a result  of  its  current  lack  of  a mechanism  for  direct  side  force 
control,  the  X-22A  aircraft  is  not  in  general  capable  of  simulating  the  lateral- 
directional  stability  and  control  characteristics  of  another  VTOL  aircraft  in^ 
toto.  However,  as  will  be  demonstrated  in  this  subsection,  the  important 
lateral-directional  flying  qualities  parameters  of  the  basic  AV-8A  may  be  very 
adequately  simulated  by  the  X-22A  using  only  the  existing  VSS  roll  and  yaw 
channels . 


The  1 ateral -directional  rigid-body  equations  of  motion  of  the  AV-8A 
may  be  described  by  a three-degree-of-freedom,  fourth-order  system  (Section  2): 

X.  - Fx  *■  G(jl 

x.  = state  vector  of  aircraft  (4  elements) 

F = aircraft  system  matrix  (3  x 4) 

G = aircraft  control  matrix  (3  x 2) 

U = aircraft  control  vector  (2  elements  - roll  and 

yaw  control) 

Thus,  twelve  state  coefficients  and  six  control  coefficients  are,  in  general, 
required  to  describe  fully  the  lateral-directional  dynamics  of  the  aircraft. 

With  two  independent  lateral-directional  controls,  roll  and  yaw,  available  in 
the  X-22A,  only  two  of  the  three  degrees  of  freedom  may  be  independently  con- 
trolled. Thus,  eight  state  coefficients-controllable  degrees  of  freedom  (2) 

X order  of  system  (4)  - and  four  control  coefficients-controllable  degrees  of 
freedom  (2)  X number  of  controls  (2)  - may  be  specified.  Therefore,  the  X-22A, 
using  only  roll  and  yaw  controls,  is  capable  of  reproducing  twelve  of  the 
eighteen  parameters  which  describe  the  lateral-directional  dynamics  of  the  basic 
AV-8A;  the  six  independent  coefficients  which  cannot  be  directly  controlled 
are  the  side  force  derivatives,  of  which  the  side-force-due-to-sideslip  ( Y ) 
has  the  dominant  effect  on  the  AV-8A's  lateral-directional  flying  qualities. 
Accordingly,  the  process  whose  end  product  is  an  acceptable  in-flight  simulation 
of  the  basic  AV-8A's  lateral-directional  characteristics  initially  involves  the 
selection  of  the  twelve  flying  qualities  parameters  to  be  specified  for  the 
simulation. 

A Calspan-developed  digital  computer  program  utilized  primarily  for 
lateral-directional  simulations  involving  the  USAF/Calspan  T-33  (Reference  28  ) 
was  used  to  calculate  the  appropriate  X-22A  VSS  gains  for  the  Harrier  lateral- 
directional  simulation.  The  input  requirements  for  this  program  include: 

1)  the  X-22A  dimensional  stability/control  derivatives  in  primed  form,  airspeed, 
and  angle  of  attack  at  the  flight  conditions  of  interest  and  2)  the  twelve 
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specified  AV-8A  stability/control  parameters  at  each  flight  condition  expressed 
in  terms  of  modal  characteristics  and  stability/control  derivatives.  Included 
in  the  program  output  for  each  flight  condition  are: 

• the  VSS  gains  required  to  achieve  the  specified  parameters 

• the  stability/control  derivatives  of  the  simulated  aircraft, 
i.e.,  X-22A  with  VSS  ir.  operation  with  calculated  gains 

• equations  of  motion  and  modal  characteristics  of  simulated 
aircraft 

• roots  and  transfer  functions  of  simulated  aircraft 

• time  histories  of  simulated  aircraft  response  to 
control  inputs. 

The  simulation  methodology,  dictated  in  part  by  the  structure  of  the 
computer  program,  was  to  specify  the  values  of  the  following  parameters  at  each 
flight  condition: 


|*/*U 

L'4  (-0) 

N'fl  (*0) 

^rP 

N'srp 

The  remaining  parameter  to  be  specified,  was  iterated  about  its  nominal 

value  in  an  effort  to  achieve  the  best  match  of  the  derivatives  L ^ and  N'  and 
the  flying  qualities  parameter  p 

Using  the  above  criteria,  the  VSS  roll  and  yaw  channel  gains  were 
calculated  as  a function  of  flight  condition,  specifically  X . These  gains 
are  presented  in  Table  4-1. 

As  discussed  above,  certain  stability  and  control  parameters  of  the 
basic  AV-8A  model  were  not  specifically  matched  by  the  VSS  gain  calculation 
techniques.  Table  4-2  presents  the  values  of  the  parameters  for  the  X-22A 
model  with  the  VSS  in  operation  at  the  calculated  gain  levels  and  compares 
them  to  the  corresponding  value  for  the  basic  AV-8A  model  described  in  Section 
II.  In  addition  to  the  differences  in  mentioned  previously,  inspection  of 
Table  4-2  reveals  that  discrepancies  between  the  characteristics  of  the  basic 
AV-8A  model  and  those  of  the  X-22A  simulation  of  the  model  exist  in  the  fol- 
lowing areas: 


TABLE  4-2 

COMPARISON  OF  STABILITY/CONTROL  PARAMETERS 
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(1)  X-22A/VSS  Simulation  of  AV-8A 

(2)  Basic  AV-8A  Model 


roll  damping  ( /.  ) 

coupling  derivatives  (Lr  and  Np ) 


A brief  evaluation  of  the  effects  of  these  discrepancies  was  conducted 
in  the  frequency  domain  by  a comparison  of  the  transfer  function  numerators  of 
the  two  models  (Table  4-3)  and  in  the  time  domain  by  a comparison  of  time 
histories  of  state  variable  responses  to  identical  roll  and  yaw  control  inputs 
(Figure  4-1).  Both  the  transfer  functions  and  the  time  history  data  indicate 
acceptable  matches  of  the  state  variables  to  pilot  control  inputs  with  the  pos- 
sible exception  of  the  yaw  rate  response  to  lateral  stick  input  in  hover.  How- 
ever the  magnitude  of  the  lateral  acceleration  response  to  control  of  the  basic 
AV-8A  model  is  significantly  less  than  the  corresponding  X-22A/VSS  response  be- 
cause of  the  order  of  magnitude  difference  in  the  Y^,  derivative  between  the  two 
vehicles.  This  same  discrepancy  will  also  result  in  widely  different  lateral 
acceleration  responses  to  turbulence. 

4.2  LONGITUDINAL 

4.2.1  Summary  of  the  Problem 

The  longitudinal  simulation  problem  is  similar  to  but  more  difficult 
than  the  lateral/directional  situation.  Again,  three  degree  of  freedom  motion 
is  to  be  simulated  with  two  controllers,  if  possible;  although  the  use  of  the 
elevons  in  a collective  fashion  would  provide  three  controllers  and  hence  exact 
simulation  could  be  possible,  the  elevons  would  have  limited  authority  and  hence 
is  was  judged  desirable  to  avoid  using  feedbacks  to  them.  Unfortunately,  the 
digital  computer  design  procedure  used  to  compute  the  lateral/directional  feed- 
back gains  was  not  an  appropriate  tool  for  the  longitudinal  design  problem 
because : 

• The  collective  blade  controller  changes  from  being  effective 
at  changing  vertical-force  equation  characteristics  at  low 
speeds  to  changing  longitudinal-force  equation  characteristics 
at  high  speeds.  Hence  a uniform  set  of  desirable  modal 
characteristics  cannot  be  simulated  over  the  entire  velocity 
range. 

• Although  the  majority  of  important  modal  characteristics  for 
lateral/directional  flying  qualities  can  be  achieved  through 
control  of  the  yaw  and  roll  moment  equations,  for  the 
longitudinal  problem  the  individual  force  equations  (particularly 
longitudinal  force)  are  not  as  effective.  In  addition, 
equivalent  modal  characteristics  (e.g.  \ $!(S  1^)  are  not  as  well 
defined  in  terms  of  flying  qualities  for  the  longitudinal 
problem,  particularly  in  transition. 
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TABLE  4-3 

LATERAL-DIRECTIONAL  CHARACTERISTIC  ROOTS  AND  TRANSFER  FUNCTIONS 


+++++++  X-22A  Simulation 

AV-8A  Model 
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Figure  4-1  LATERAL-DIRECTIONAL  TIME  HISTORIES  (a) 


Figure  4-1  LATERAL-DIRECTIONAL  TIME  HISTORIES  (f)  V = 120  kt;<£rp  DOUBLET 


( 
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As  a result  of  these  difficulties,  alternate  design  procedures  were 
investigated  for  the  longitudinal  gain  calculations.  In  particular,  "least- 
squares"  and  optimal  control  techniques  were  examined  for  applicability  to  the 
feedback  gain  computations,  while  least  squares  and  modified  response-error 
techniques  were  investigated  for  the  gearing  calculations.  The  following  sub- 
section presents  a tutorial  review  of  the  feedback  design  procedures,  indicates 
the  reasons  for  selecting  the  optimal  control  procedure,  and  develops  a 
theoretical  extension  to  the  optimal  control  methodology  that  is  required  for 
this  application.  The  feedback  design  results  are  then  presented  and  discussed 
in  Section  4.2.3.  Possible  ways  of  computing  the  gearings  are  discussed  and 
developed  in  Section  4.2.4  along  with  the  results  obtained  with  the  selected 
method.  Finally,  time  history  responses  of  the  resulting  longitudinal  simula- 
tion are  compared  with  the  AV-8A  model  responses  in  Section  4.2.5. 

4.2.2  Development  of  Feedback  Laws 

The  least-squares  result  is  obtained  by  minimizing  a non-integral 
performance  index  of  the  general  form: 

V = erQ  e *■  urR  a. 

This  problem,  it  should  be  observed,  is  a direct  minimization  procedure  since 
no  integral  is  involved:  to  find  the  value  of  the  control  which  minimizes  V 

we  set  the  gradient  (derivative)  of  V with  respect  to  u equal  to  zero: 


V V - O , solve  for  Li. 
u ’ 

Since  the  intent  of  the  design  is  to  compute  feedback  gains  which  will  make 
one  system  (the  "plant")  match  the  characteristics  of  another  system  (the 
"model") , a convenient  choice  for  the  error  term  in  the  performance  index  is 
the  equation  error  between  the  two  systems.  To  concentrate  our  attention  on 
the  feedback  computations,  consider  the  following  example: 

x - ■+  G a (Plant) 

■ ij  ^ Hz 

y r L y (Model) 

so  that  e = z - y = { f - L H)  x + G u. 

Choosing  H-l  for  simplicity  (i.e.,  the  model  states  correspond  directly  to 
the  plant  states),  and  substituting  into  the  equation  for  ]/  , we  find: 

1/  = xr(F~L)  T<?  (P-£  )x  *-  2xr(F-L)rQCu  + ur  [orQG  + Z]u. 

Then: 

ZGtQ  (F-L)x  + 2 \GTQG  * *]oi 
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The  optimal  control  result  is  obtained  by  minimizing  an  integral  per- 
formance index  of  the  general  form: 

I = f (crOe  + oc T R u. ) cLt 
Jo 

The  problem  now  is  the  minimization  of  a functional,  which  requires  the  use  of 
somewhat  more  complicated  mathematics,  such  as  the  calculus  of  variations  or 
Pontryagin's  maximum  principle,  to  solve.  The  problem  has  been  solved  and 
discussed  for  many  choices  of  the  error  term  (see  References  29,  30);  although 
repeating  the  development  here  is  not  therefore  strictly  necessary,  it  may  be 
a useful  exercise  to  understand  the  result  obtained.  Accordingly,  a brief 
development  of  this  particular  example  follows,  with  some  liberties  taken  in 
the  cause  of  simplification. 

For  the  same  error  expression  used  in  the  least-squares  example, 
our  problem  is: 

Minimize  *00  r 

X =■  J [_(z  - y)T Q(x  - y ) + ct  Ru]dt 

0 

- J [(Fx  + Gu.  - Lx)1 Q(Fx  + Gu-  Lx)  + U.T(?u.'\ctt 


subject  to  * = Fx  + G u. 

Recall  that,  in  the  calculus  of  variations  (to  use  the  solution  method  most 
familiar) , the  system  constraints  are  introduced  through  the  use  of  Lagrange 
multipliers : 

/ V Q J'  [(Fx  + Gu.  + Lx)7 Q ( Fx  + Gu  + Lx)  + u.T Ru  + 2 A • (~z  + Fx  + Gu 

A necessary  condition  for  the  optimum  to  exist  is  then  that  the  Euler  equations 
plus  the  constraint  equation  be  satisfied.  These  relationships  are: 
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ci  r jj  \ ££ 

ait  \ d LL  ) d U.  ^ 

X -Fz  - Gu  - O 

Considering  the  second  equation: 


where  J is  the  integrand 


14-0 

3 OL 


(u.  does  not  appear  explicitly) 


Hence : 


~ ^ ^ rQ  ( F x.  + G u.  - lx)  ■+•  2 Rcjl  2 G 7" 2. 


2GrO  (F-L)x  + 2(GrOG  + R)  u.  + 2 GrA=0 


or: 

= - [gTQG  + [GTQ  (F-L)x  + GT *] 

Before  continuing  with  the  solution,  it  can  be  immediately  noted  that  the  form 
of  the  result  for  the  control  is  different  than  the  one  obtained  with  the  least- 
square  computation:  the  Lagrange  multiplier  (or  "adjoint  variable")  is  added 

to  the  least-square  result. 


Hence : 


Considering  now  the  first  equation: 


3T 
3 x 


-2  A , 


dJ 


2{F  -L)T Q (Fx  -r  Gu.  - Lx) 


2Fr-\ 


- \ - FrX 


(F  - L)TG  (Fx  + Gu  - Lx)  = 


We  now  have  three  equations  in  the  three  unknowns  ( x , a. , "X  ) . Eliminating 
ix  by  using  the  solution  of  the  second  equation  we  have: 


"A  - £ F - ( F - L)  Q Cr  QGr  + r\  (j  j A 

-{( F~L  )TQ(F-L)  - (F- L)TQG  [g TQG+  r]gTQ(F-l)\x  = O 

x - \f  - G[GrQG  + r]  1GTQ(F'L)\x  * & \_Gr  QG  + Gr  * = 0 


These  two  equations  constitute  the  necessary  conditions  for  the 
optimum  control  to  exist.  It  is  necessary  to  solve  them  for  A in  order  to 
obtain  the  optimum  control,  as  is  evident  in  the  expression  for  aOPr  . To 
obtain  this  solution,  it  is  easiest  to  define  first  some  expressions  to  simplify 
the  form  of  the  equations;  although  these  definitions  seem  arbitrary  at  first 
glance,  the  point  is  to  put  the  equations  in  the  "canonical"  form  of  the  regu- 
lator problem  ( e = x in  the  performance  index),  for  which  the  solution  is  known. 


Let : 


F = F - G \_(jTQ  G + R~\~7  GTQ  ( F - L) 


Q = Q - QG\GrQG  A R]~'g  TQ 


R t 6TQCx 


H = F -L 


Then  the  equations  are: 


a r a * 


" A - F A +■  AV  Q A/  x 


. A ^ < f 

x = F x - GR  6 a 


We  now  assume  that  the  adjoint  variable  has  the  form: 
A = P(t)x 

Then: 


Px  + Px.  = Px  + PF  x - PGR~f  QrPx 
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Since 


X 


Fx  * 


A - / r 
GR  G A 


A A - / T 

Fx  - GR  G Px 


Substituting  into  the  A equation: 

Px  + PFx  - PGR~’  G Px  + F Px  * F QFx  = 0 

The  assumption  on  the  form  of  ^ means  that  this  expression  must  hold  for  all 
x ; hence: 

- A A *-  A . ~r  S'-  7 A T 

P + PF  + F P - PGR  ~ G P = Q H 


This  equation  is  called  the  Ricatti  equation,  and  is  a nonlinear  differential 
equation  in  P . The  use  of  infinity  as  the  upper  integration  limit  in  the  per- 
formance index,  however,  permits  us  to  consider  only  the  steady-state  solution, 
which  is  found  from  the  algebraic  equation  obtained  by  setting  P - 0.  Finally, 
then,  we  have  the  optimal  control  solution  as: 

^opt  = -\GrqG  *r]~,gt  \q(f-l)+p\x 


where 

A / T ■F'  /\  /'A 

PP  F P - PGR-  G p = -H  QH 

Comparing  the  optimal  control  and  least-squares  solutions,  the  dif- 
ference in  the  designed  feedback  gains  is  clear.  This  difference  is  related  to 
the  following  fundamental  characteristics  of  the  two  procedures: 

• The  optimal  control  procedure  guarantees  a stable  closed-loop 
system*  regardless  of  whether  or  not  the  plant  (or  model)  is 
unstable  and  regardless  of  the  number  of  controllers  used.  No 
such  guarantee  exists  for  the  least-squares  solution. 

• The  least-squares  result  will  reproduce  the  model  characteristic 
roots  only  if  three  controllers  are  used.  For  one  or  two 
controllers,  some  of:  the  closed-loop  system  roots  may  be 
unstable  even  if  the  model  is  stable.  If  the  model  is  stable, 
the  optimal  control  result  can  provide  characteristic  roots 
close  to  those  of  the  model  (depending  on  the  relative  Q and 

R weightings)  for  less  than  three  controllers. 

The  desirability  of  using  only  two  controllers  makes  the  least- 
squares  technique  unacceptable  for  this  application.  The  optimal  control  pro- 
cedure is  inherently  more  suitable  to  situations  in  which  complete  controllability 
(3  controllers)  does  not  exist.  However,  another  characteristic  of  the  optimal 
control  causes  a different  problem  for  this  application:  the  closed-loop  is 


* Subject  to  R being  positive  definite,  Q positive  semi-definite 
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guaranteed  stable,  but  the  AV-8A  model  is  unstable.  Hence,  the  roots  of  the 
augmented  plant  using  the  optimal  control  feedback  gains  do  not  include  the 
unstable  model  roots,  and  would  therefore  provide  a poor  simulation.  It  was 
therefore  necessary  to  extend  the  plant-in-the-performance-index  optimal  con- 
trol theory  to  deal  with  this  situation. 

As  has  been  discussed  by  Rynaski  (Reference  29),  a useful  visualiza- 
tion of  the  optimal  control  procedure  is  the  root-square  locus  technique.  This 
technique  utilizes  the  characteristics  of  the  Hamiltonian  to  show  that  the 
adjoint  roots  are  just  the  reflection  of  the  system  roots  about  the  imaginary 
axis,  and  that,  with  both  sets  of  roots  drawn,  conventional  root-locus  tech- 
niques can  be  used  to  compute  their  migrations.  When  the  model-in-the- 
performance- index  procedure  is  used,  the  roots  of  the  plant  are  the  poles  and 
the  roots  of  the  model  are  the  zeroes.  The  loci  in  the  left-hand  plane  are 
then  the  optimal  control  solutions  for  varying  <?  and  R . 

Two  examples  of  root-square  locus  sketches  are  shown  in  Figure  4-2 
for  a second-order  system:  in  (a),  both  the  plant  and  model  are  stable, 

whereas  in  (b)  the  plant  is  stable  but  the  model  is  unstable.  The  adjoint  poles 
and  zeroes  are  shown  as  dashed  symbols. 


\ 

i JW 

l 

c 

? 

V 

1 

1 

% 

/ 

? 

1 

“^^***IB 

1 

1 

'x  1 r 

1 cr 

1 

1 <T" 

| 

l 

) 

i 

JL 

' * *> 

• 

1 

6 

Optimal 

! Optimal 

Control 

Control 

Solution 

Solution 

a)  Stable  Model  b)  Unstable  Model 


Figure  4-2  OPTIMAL  CONTROL  EXAMPLES  IN  S-PLANE 
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As  can  be  seen  from  these  sketches,  what  happens  when  the  model  is  unstable  is 
that  the  plant  roots  migrate  to  the  adjoint  zeroes  (which  are  stable)  to  ensure 
a stable  augmented  system;  in  fact,  for  this  example  the  optimal  gains  would 
be  computed  to  be  the  same  for  both  (a)  and  (b) . The  optimal  control  procedure 
will  always  use  the  poles  and  zeroes  on  the  left-hand  side  in  computing  the 
gains,  and  will  never  try  to  "follow"  an  unstable  model  if  Q and  R meet 
positive-definite  requirements. 

To  use  the  optimal  control  procedure  for  our  application,  therefore, 
it  has  been  necessary  to  devise  an  approach  to  "fool"  the  computational  tech- 
nique. The  oasic  idea  is  to  perform  the  analysis  in  a transformed  s-plane 
which  is  selected  such  that  both  plant  and  model  are  stable  with  respect  to  it; 
the  gains  which  give  desired  transformed  eigenvalues  are  then  shown  to  yield 
following  of  the  unstable  model  in  the  original  s-plane. 

The  eigenvalues  of  the  plant  and  model  can  be  expressed  by: 

P'f FP  = Ap,  n>-'L  M --  Am 

■ ^ 

Choose  a positive  number  "n"  such  that,  when  it  is  subtracted  from  the  most 
unstable  eigenvalue,  real  or  complex,  the  new  eigenvalue  is  stable.  Then  define: 

^ p ~ V ~ ^ I y m ~ ~ ^ ^ 

It  is  then  easy  to  show  that: 

r-'f'p  r if  r'  = F-nl 

M ''l'hi  - Al m if  L - L - ml 

We  now  solve  the  following  optimal  control  problem: 


Minimize 

j-  - - y)r  q (*  - y)  + 

Subject  to:  • _/  , 

X = F x + F>u. 

y = L y 
* = y 

For  this  transformed  problem,  all  of  the  system  poles  and  zeroes  are  on  the 
left-hand  side  of  the  transformed  s-plane  and  all  the  adjoint  roots  are  in  the 
right-hand  side;  hence,  conventional  juggling  of  the  Q and  R matrices  can 
be  used  to  get  the  augmented  plant  roots  (in  the  transformed  plane)  as  close 
to  the  model  roots  as-  desired. 
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Calling  the  resulting  optimal  control  gain  matrix  K and  the  augmented 
eigenvalues  in  the  transformed  system  , we  may  now  w rite: 

A ( F - & K ) A ~ JV.  A 

Then: 

A'1  (F  - GK')  A = JLA  + nl 


Hence,  using  the  computed  gains  with  the  original,  untransformed  plant  matrix 
yields  a set  of  eigenvalues  A4  1 if  the  optimal  control  procedure  made 

-A  4 = A-m*  t^len  A- A ~ A m as  desired  because  V\.m  = + nl'  Figure  4-3 

sketches  the  result  of  this  procedure  for  the  unstable  model  example  shown  in 
Figure  4-2b. 


Transformed  S-Plane 

Figure  4-3  EXAMPLE  OF  TRANSFORMED  OPTIMAL  CONTROL  SOLUTION 


With  this  modified  technique,  it  is  possible  to  use  optimal  control 
procedures  to  compute  feedback  gains  to  simulate  the  AV-8A  model  in  the  X-22A. 
The  following  subsection  presents  the  results  of  these  analyses. 
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4.2.3  Feedback  Gain  Computations 

Gains  for  the  longitudinal  simulation  of  the  AV-8A  were  computed  at 
the  six  flight  conditions  discussed  in  Section  II:  V0  - 0,  30,  60,  80,  100, 

120  kts;  y - -S  degrees,  cc  = +8  degrees.  The  model  values  used  in  the 
model  matrix  were  given  in  Table  2-1;  the  plant  values  are  the  X-22A  stability 
and  control  derivatives  programmed  on  the  ground  simulator  as  listed  in 
Table  5-1  in  the  next  section.  For  the  X-22A,  trim  velocity  and  duct  angle 
are  closely  related;  hence,  the  computed  gains  may  be  programmed  in  the 
variable  stability  system  as  a function  of  either  variable.  For  this  applica- 
tion, duct  angle  was  selected  since  the  signal  is  monotonic  and  less  noisy 
than  the  velocity  signal:  the  computed  gains  to  be  discussed  here  will  there- 

fore be  presented  for  flight  condition  in  terms  of  duct  angle  as  well  as  trim 
velocity. 


The  gains  for  each  flight  condition  were  computed  assuming  two  con- 
trollers: longitudinal  stick  and  collective  blade  pitch.  Since  only  two  con- 

trollers are  used,  perfect  matching  of  all  the  model  characteristics  is  impos- 
sible, and  hence  it  is  necessary  to  decide  which  characteristics  are  most  im- 
portant to  match.  For  this  application,  matching  the  characteristic  roots  of 
the  model  matrix  was  judged  most  important  --  in  particular,  matching  the  most 
stable  and  most  unstable  roots  was  of  paramount  importance.  The  values  of 
feedback  gains  selected  were  the  ones  which  best  achieved  this  desideratum. 

The  computational  procedure  consists  of  selecting  various  values  of 
the  O and  £ matrices,  computing  the  gains,  and  checking  the  resulting  closed- 
loop  eigenvalues.  This  iterative  "art"  was  a fairly  straightforward  process 
for  the  low  speed  configurations,  because  the  two  controllers  have  major 
effectivenesses  on  the  moment  and  vertical  force  equations;  the  higher  speed 
cases  ( V0  = 80,  100,  120  kt)  were  more  difficult  because  the  collective  pitch 
effectiveness  is  primarily  on  the  longitudinal  force  equation.  For  these  latter 
cases,  it  was  necessary  to  assume  different  values  of  Q and/or  R in  the 
model  matrix  to  achieve  feedback  gains  which  made  the  closed-loop  roots  as 
unstable  as  the  original  model  roots.  A summary  of  the  values  of  O and  R used 
is  given  in  Table  4-4. 

The  computed  feedback  gains  are  given  in  Table  4-5a,  and  simplified 
(for  implementation  purposes)  values  listed  in  Table  4-5b.  The  closed-loop 
roots  using  the  simplified  gains  are  compared  with  the  model  roots  in  Table 
4-6.  As  can  be  seen,  the  correspondence  is  generally  good  for  the  dominant 
roots.  Table  4-7  compares  the  closed-loop  characteristic  matrices  with  the 
model  matrices;  note  that  the  M and  B equations  are  well  matched  for  VQ  = 0, 

30,  60  kts,  but  that  the  X and  £ equations  are  dissimilar  for  \XQ  = 80,  100, 

120  kts.  This  dissimilarity  leads  o poor  matching  of  time  history  responses 
to  throttle  inputs,  as  will  be  discussed  shortly. 
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TABLE  4-4 


WEIGHTING  MATRICES  AND  MODEL  MODIFICATIONS 
FOR  LONGITUDINAL  GAIN  COMPUTATIONS 


V 

. «i 

R 

Modified  Model 

0 

diag  (.001,  10,  .1,  100) 

diag  (.3,  .1) 

None 

30 

diag  (.001,  10,  .1,  100) 

diag  (.3,  .1) 

None 

60 

diag  (.001,  10,  .1,  100) 

diag  (0,  0) 

None 

80 

diag  (.001, .001, .001,100) 

diag  (.3,  .1) 

= .0051,  = -.953 

100 

diag  (.001, .001, .001, 100) 

diag  (.3,  .1) 

= . 00545 , Mp  = -1.06 

120 

diag  (.001, .001, .001, 100) 

diag  (.3,  .1) 

A*,=  .006,  = -1.35 

TABLE  4-5a 

LONGITUDINAL  FEEDBACK  GAINS: 
EXACT  COMPUTATIONS  IN  COMPUTER  UNITS 


D 

As/u 

Ae/U r 

Ae/G 

Ae/9 

Ac/e 

Ac/<Z 

0 

.039 

0 

.0007 

+ .633 

-.0012 

+ .066 

0 

-.0672 

30 

-.028 

-.044 

-.109 

+ .603 

+ .113 

+ .271 

2.56 

-.0491 

60 

-.053 

-.101 

- .064 

+ .373 

+ .197 

+ .419 

+ .841 

+ .020 

80 

-.032 

-.049 

-.016 

- .964 

-.012 

-.015 

.0084 

-.44 

100 

-.028 

-.032 

-.013 

-1.25 

-.011 

-.0089 

- .0053 

-.525 

120 

-.026 

-.031 

-.011 

-1.02 

-.013 

-.010 

- .013 

-.510 
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TABLL  4-5b 


LONGITUDINAL  GAINS  FOR  IMP  LIMITATION 
IN  GROUND  SIMULATOR 


^fS  /U  AfS/ur  Azs/O 


-.0111 
-.0105 
- .0065 
+.0168 
.0147 
.0178 


^ClA^  ^ CS /“r  Acs  /Q  Acs/<2 


- . 06b  0 

-.271  -.0440 

-.419  -.0147 

+ .012  +.015  0 

+ .011  0 0 


0 

0 

0 

.0077 

.0066 

.0089 


TABLL  4-6 

COMPARISON  OF  MODLL  AND  SIMULATION 
CHARACTERISTIC  ROOTS  IN  FORM  ( l/Z) [ X i 6d„  J 


mm 

AV-8A  Model 

X-22A  Simulation  of  AV-8A 

0 

( . 33)  ( . 02  ) [-.48;  .3lJ 

( .40)  (.02)  £.40;  .33] 

30 

( .66) (.064)  [-.70;  .33  J 

( .61)  (.10)  [-.42;  .28] 

60 

(1.00)  (.073)  [-.90;  .32  ] 

(1 . 05) (.15)  [-.83;  .25] 

80 

(1 .21) (.075)  (- . 38) (- . 26) 

(1.14) (.16)  (-.37) (-.12) 

100 

(1 . 39) ( . 079)  (- . 48) (- . 20) 

(1 .37)  (. 16)  (-.43)  (-.062) 

120 

(1 .53) (• 084)  (- . 47) (- . 19) 

(1.46) (. 12)  (- . 47)  ( .078) 
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TABLE  4-7 


COMPARISON  OF  MODEL  AND  SIMULATED 
CHARACTERISTIC  MATRICES 


V 

AV-8A  Model 

X-22A  Simulation  of  AV-8A 

0 

-.026  0 -32.2  0 

0 -.020  0 0 
0 0 0 1.0 

.001  .0009  0 0 

-.145  0 -32.2  -2.91 

.0018  -.020  0 .10 

0 0 0 1.0 

.00135  .00092  0 .0045 

30 

-.035  -.020  -32.1  -7.1 

-.011  -.105  - 1.68  50.7 

0 0 0 1.0 
_ .0017  .0042  0 - .13 

-.222  -.156  -33.07  -11.81 

-.012  -.107  - 1.69  50.75 

0 0 0 1.0 
_ .00129  .0030  - .035  - .154  j 

60 

-.043  -.027  -32.1  -14.2 

-.036  -.19  - 1.68  101.0 

0 0 0 1.0 
.0022  .0054  0 - .26 

U.  - 

-.333  -.34  -32.67  -18.87 

-.036  -.19  - 1.68  101.0 

0 0 0 1.0 
.0022  .00526  - .020  - .26 

80 

1 

BHDH 

H 

100 

-.055  -.030  -32.1  -23.8 

-.083  -.31  - 1.68  170.0 

0 0 0 1.0 

.0026  .006  0 - .43 

L _/ 

-.173  .092  -32.1  -29.45 

-.258  -.639  - 2.81  170.3 

0 0 0 1.0 
.0022  .0052  0 - .231 

120 
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4.2.4 


Control  Effectiveness  Computations 


A major  difference  between  the  AV-8A  model  and  the  X-22A  is  the  thrust 
incl  nation  as  a function  of  velocity,  both  for  trim  conditions  and  particularly 
for  the  deceleration  profile.  This  difference  implies  difficulty  in  simulating 
ver  ical  and  longitudinal  acceleration  and/or  velocity  responses  to  throttle 
inputs  at  higher  speeds  if  only  collective  blade  pitch  and  longitudinal  stick 
controllers  are  used  in  the  X-22A.  For  this  reason,  the  use  of  the  X-22A 
elevons  in  a collective  fashion,  which  requires  a minor  modification  to  the 
X-22A  variable  stability  system,  appears  to  be  required  for  the  simulation  of 
the  AV-8A  control  effectiveness  characteristics;  this  capability  has  been 
assumed  in  the  computations. 

Several  methods  for  computing  the  control  effectiveness  gearings  were 
investigated.  The  two  that  offered  the  most  promise  and  received  the  most  atten- 
tion were  least-squares  equation  error  and  weighted  response-error  minimization. 
The  former  is  the  control  matrix  equivalent  of  the  least-squares  feedback  compu- 
tational procedure;  the  latter  is  an  extension  of  the  method  used  for  the  design 
of  the  decoupled  velocity  control  system  in  the  Task  III  X-22A  experiment 
(Reference  32) . 

Th  least  squares  result  is  obtained  by  minimizing: 

J * e.rQe  , e = Gu  - Um 


Then : 


U = (GrQG  )~fG  rQ  Gm  U.m 

This  method  will  result  in  the  simulated  control  derivatives  being  exactly  the 
same  as  the  model  derivatives.  For  all  situations,  therefore,  the  initial  ac- 
celerations are  properly  matched.  If,  in  addition,  the  characteristic  matrix 
has  been  exactly  matched,  then  a one-for-one  simulation  results;  as  can  be 
seen,  however,  the  simulated  characteristic  matrix  does  not  enter  into  this 
computation,  and  if  it  is  not  well  matched  (such  as  for  the  VQ  = 80,  100,  120 
kt  cases),  the  time  history  responses  after  zero  time  may  not  be  a good  simula- 
tion of  the  model  responses. 

The  method  used  to  design  the  decoupled  velocity  control  system  gear- 
ings was  an  attempt  to  assure  correct  time  history  responses  (in  terms  of  steady 
state  magnitudes)  even  if  the  characteristic  matrix  is  not  totally  matched.  The 
closed-loop  equation  is  written  as: 

X -*  F.  X + GJ  um 
where  J~  are  the  control  gearings. 

For  unit  inputs,  the  steady  state  responses,  assuming  a stable  system,  are: 
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I 


X 


= - F-'GJ  = - PJ 


Minimizing  the  difference  between  these  responses  and  the  desired  steady-state 
responses  (X^,)  leads  to  the  equation  for  the  gearings: 

J --  - (PrQ  pY'  PrQ 

Although  this  procedure  was  found  to  be  the  most  suitable  for  the 
Task  III  design,  it  was  not  directly  applicable  to  the  AV-8A  simulation  problem 
because  the  closed-loop  system  is  not  stable.  Accordingly,  steady-state 
responses  do  not  exist  and  the  solution  X = - PJ  is  not  valid.  An  attempt  was 
therefore  made  to  extend  the  analysis  to  unstable  systems  as  well. 

The  state  responses  at  any  time  to  unit  step  inputs  are,  regardless 
of  whether  the  system  is  stable  or  unstable: 

x = re~'(x  - oj) 

Then: 

XD-X  = x0  - Fc-'(x  - GJ) 

• < 

If  we  assume,  in  order  to  minimize  this  error  term,  that  we  desire  X = XQ  also, 
then  the  expression  for  J"  becomes: 

s - - (pt«p)''  pTQ[xD’FcXil 

This  expression  was  used  to  investigate  gearing  computations  at  the  60  kt  flight 
condition.  As  can  be  seen  from  examining  the  equation,  a major  difficulty  is 
that  the  results  for  the  gearings  will  be  different  for  each  time  selected  in 
obtaining  Xp  , X^j  ; hence,  unless  a criterion  exists  which  specifies  that  the 
responses  shall  be  the  same  at  one  particular  time,  no  rational  method  of  se- 
lecting the  "best"  gearings  on  this  basis  is  apparent. 

The  least-squares  equation-error  procedure  was  therefore  used  to  com- 
pute the  gearing  gains;  the  results  are  listed  in  Table  4-8.  Since  three  con- 
trollers are  assumed,  these  gains  reproduce  the  AV-8A  model  control  matrix 
exactly.  The  initial  accelerations  will  therefore  be  correct,  but  no  guarantee 
about  the  remainder  of  the  time  history  responses  can  be  made:  differences 

will  be  dependent  on  differences  in  the  characteristic  matrices. 
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TABLE  4-8 

LONGITUDINAL  GLARING  GAINS 


I 


v/ 

&es/  &T 

^ cs  / &T 

■SKA 

&6L  /Se 

A £ i / &t 

mmmammm 

0 

B 

-.0018 

-.32 

2.12 

- .15 

- .95 

30 

B 

-.087 

-.37 

2.50 

- .93 

1.66 

60 

.80 

-.135 

-.22 

2.22 

-1.54 

4.20 

80 

.91 

-.147 

-.023 

1.64 

-1.51 

5.70 

100 

1.06 

-.060 

-.174 

.95 

- .38 

6.84 

120 

1.20 

.023 

-.178 

.28 

- .26 

7.78 

00 

O 

* 

.91 

-.178 

-.023 

2.04 

-1.51 

4.49 

120** 

1.20 

-.160 

-.178 

1.48 

- .26 

4.91 

* 03 

o 

o 

r- 

ii 

No  asterisk: 

3 - «° 

**  0J 

= 50° 
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4.2.5  Resulting  Longitudinal  Simulation 

Time  history  responses  to  simulated  longitudinal  stick  and  throttle 
pulse  inputs  are  given  in  Figure  4-4  for  the  V0  = 0,  60,  120  kt  trim  conditions. 
The  solid  lines  are  the  AV-8A  model  responses,  and  the  crosses  are  those  of  the 
X-22A  simulation  of  this  model.  As  can  be  seen,  excellent  matching  of  all  state 
variable  responses  to  pitch  stick  inputs  is  achieved  for  the  low  speed  (0,  60  kt) 
cases;  the  30  kt  responses  (not  shown)  are  similar.  The  only  real  discrepancy 
is  in  the  nx  reponses,  which  results  fro  the  dissimilar  values  of  xu  and  Xv\ 
these  derivatives  cannot  be  modified  significantly  at  the  low  speeds  (high  duct 
angles)  with  the  collective  plade  pitch  controller,  and  matching  nK  exactly  would 
require  a third  controller.  Hence,  for  example,  the  longitudinal  velocity  re- 
sponse (a)  to  throttle  inputs  at  60  kt  (Figure  4-4d)  is  not  a good  match,  but 
the  magnitudes  are  still  small  relative  to  the  well-matched  vertical  velocity 
responses.  It  is  not  considered  that  this  discrepancy  will  have  a significant 
flying  qualities  effect  on  the  simulation. 

Good  matching  of  the  state  variable  responses  to  pitch  stick  inputs 
is  also  apparent  at  the  120  kt  case;  similar  matches  were  obtained  at  80,  100 
kt  (not  shown).  Discrepancies  exist  for  both  nx  and  because  neither  X - 
force  nor  2 -force  derivatives  are  replicated  at  the  higher  speeds;  the  inability 
to  vary  £u  and  with  two  controllers  when  the  duct  angle  is  low  is  reflected 
in  the  differences.  At  120  kt,  good  matches  of  the  responses  for  throttle 
inputs  were  obtained  also  (Figure  4-4f),  but  the  vertical  velocity  responses  at 
80  and  100  kt  (not  shown)  were  considerably  different.  This  difference  is  a di- 
rect result  of  the  inability  to  alter  the  X-22A's  2 ^ derivative  at  low  duct 
angles  without  the  use  of  a third  controller,  coupled  with  the  selection  X(,T 
and  Z&T  to  be  the  same  as  for  the  AV-8A  model.  The  effect  of  these  mismatches 
can  be  seen  by  considering  an  approximate  version  of  the  transfer  function 
numerator  of  vertical  velocity  to  throttle: 

"Z  * ts,  s 


q * * 


for  S*&r~  O,  Qo  = o. 

The  situation  is  analogous  to  the  placement  of  the  zeroes  in  the  roll-to-lateral- 
stick  transfer  function  using  yaw-due-to-ai leron  ( /Vj^) : The  Xsr  /is7  ratio 

may  be  changed  to  alter  the  zero  locations  in  this  transfer  function  at  the 
expense  of: 

(1)  Non-matched  initial  acceleration  responses 

(2)  Variations  in  the  longitudinal-velocity-to-throttle  transfer 
function  numerator. 
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It  is  not  clear,  however,  the  extent  to  which  these  zeroes  should  be  matched  if 
the  longitudinal  velocity  transfer  function  zeroes  become  less  well  matched 
simultaneously.  This  discrepancy  will  therefore  require  additional  attention 
during  the  initial  stages  of  the  ground  simulation  experiment  in  the  follow-on 
flight  program. 

4.3  SUMMARY  OF  SIMULATION  FIDELITY 

This  section  has  discussed  the  procedures  used  to  design  the  gains 
required  to  simulate  the  AV-8A  model  with  the  variable  stability  X-22A.  From 
0 kt  to  60  kt,  good  duplication  of  the  state  variable  responses  for  all  control 
inputs  was  generally  achieved,  with  some  reduction  in  longitudinal-velocity-to- 
throttle  fidelity  at  60  kt.  The  only  discrepancies  for  these  flight  conditions 
are  in  the  longitudinal  and  lateral  acceleration  reponses.  Comparison  of  flying 
qualities  parameters  from  MIL-F-83300  for  the  AV-8A  model  and  the  simulation  pre- 
dict equivalent  flying  qualities.  From  80  kt  to  120  kt,  the  lateral-directional 
duplication  remains  excellent  (except  for  n y ),  and  duplication  of  longitudinal 
responses  to  pitch  stick  inputs  is  good  for  the  state  variables;  neither  the 
longitudinal  nor  vertical  acceleration  responses  match  perfectly.  Vertical 
velocity  responses  to  throttle  inputs  are  considerably  different,  however,  at 
80  kt  and  100  kt;  it  is  anticipated  that,  at  120  kt  at  least,  the  throttle 
will  be  primarily  a speed  controller  and  so  this  discrepancy  will  not  be  a 
major  deficiency.  It  is  recommended  nonetheless  that  further  analyses  of  this 
discrepancy  be  conducted  and  directed  at: 

• Alternate  simulation  designs:  either  limited  use  of  a 

third  controller  in  the  feedbacks,  or  variation  of  /Z& 
away  from  the  AV-8A  model  value.  r 

• Pilot-in-the-loop  analysis:  attempt  to  ascertain  likely 

control  techniques,  assess  impact  of  discrepancy. 

• Piloted  ground  simulation:  compare  AV-8A  model  and  the 

simulation  with  piloted  evaluations. 

The  calculations  have  shown,  therefore,  that  a generally  good  simula- 
tion of  the  AV-8A  can  be  designed  for  implementation  on  the  X-22A.  The  follow- 
ing section  discusses  the  equally  difficult  task  of  implementing  this  design. 
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Section  5 


IMPLEMENTATION  IN  GROUND  SIMULATOR 


5.1  SIMULATOR  PREPARATION 

In  an  effort  to  make  the  operation  of  the  X-22A  ground  simulator  more 
like  that  of  the  aircraft  and  thus  to  produce  a more  useful  research  tool  for 
both  ground-  and  in-flight  simulation  programs,  several  modifications  were  made 
to  the  simulator  facility  as  part  of  this  preliminary  program.  They  include: 

• the  updating  of  the  basic  X-22A  mathematical  model  as 
programmed  on  the  simulator  digital  computer 

• the  integration  of  the  airborne  analog  computer  into 
the  simulator  facility 

• the  incorporation  of  the  simulator  VSS  into  the  operation 
of  the  facility. 

Based  upon  the  identification  results  of  the  X-22A  Task  III  program, 
documented  in  Reference  31,  the  original  values  of  the  stability  and  control 
derivatives  which  provide  the  basis  for  the  simulator  mathematical  model  of  the 
basic  X-22A  in  transition  were  altered  to  provide  an  updated  version  of  that 
model.  The  current  values  of  the  X-22A  stability  and  control  derivatives  are 
presented  in  Table  5-1  at  the  same  six  flight  conditions  for  which  the  AV-8A 
model  presented  in  Section  2 is  developed. 

During  the  Task  III  flight  program,  the  airborne  analog  computer  per- 
formed the  following  functions: 

• guidance  computations 

• synthesis  of  control  director  laws  and  other  display  logic 

• implementation  of  stability/control  augmentation  systems. 

A portion  of  the  preliminary  program  documented  herein  was  devoted  to  the 
incorporation  of  the  analog  computer  into  the  ground  simulator  facility  to 
accomplish  these  same  functions,  originally  performed  by  three  F.AI-380  computers. 
The  current  analog  computer  program  for  the  simulator  is  similar  to  that 
documented  in  Reference  31. 

The  requirement  to  sim-’late  a specific  VTOL  aircraft  in  transition 
for  the  upcoming  program  made  mandatory  the  use  of  the  simulator  VSS  with  its 
capability  of  feedback  gain  variation  as  a function  of  airspeed  or  duct  angle. 
Hence  the  system  was  made  fully  operational  during  the  course  of  the  preliminary 
program. 
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TABLE  5-1 

BASIC  X-22A  SIMULATOR  MODEL 


Deriv- 

VCA) 
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30(65) 

60  (45) 

80(30) 

100(15) 

120(0) 

Kts(Deg) 

X 

UL 

-.15 
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-.19 

-.19 
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.015 
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1.35 

2 
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1 

Lri 

O 

2 
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.310 
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2 
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; (a)  Longitudinal 
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The  performance  of  the  revised  digital  computer  program,  analog  com- 
puter, and  VSS  was  evaluated  by  reproducing  the  final  Task  III  evaluation  task, 
angular  augmentation  control  systems,  and  display  logic  on  the  simulator  using 
the  new  equipment  and  then  conducting  piloted  decelerating  approaches  in  order 
to  reveal  any  problem  areas  in  the  simulation. 

The  following  subsection  discusses  the  functions  of  the  analog  com- 
puter and  VSS  and  their  integration  for  the  simulation  of  the  AV-8A  in  transits 

5.2  AV-8A  MODEL  MECHANIZATION 

The  simulation  of  the  AV-8A  stability  and  control  characteristics 
developed  in  Section  IV  was  implemented  on  the  ground  simulator  in  a manner 
which  was  judged  to  be  both  logical  from  an  experimental  point  of  view  and  com- 
patible with  the  actual  aircraft  experimental  systems.  Specifically,  the  simu- 
lator digital  computer  retains  the  function  of  calculating  the  basic  X-22A 
equations  of  motion  and  the  guidance  parameters  that  will  be  supplied  in  the 
flight  program  by  the  microwave  landing  system  (MLS).  The  VSS  is  used  to  the 
extent  possible  to  provide  the  characteristics  of  the  basic  AV-8A  model;  the 
analog  computer  will  be  used  to  generate  the  models  of  the  SCAS  for  the  up- 
coming AV-8B  simulation  program  and  is  currently  programmed  with  a simplified 
model  of  the  AV-8A  SAS  (documented  in  Appendix  II).  Although  it  is  possible  to 
mechanize  many  of  the  control  system  schemes  directly  on  the  VSS  and  thus  free 
a good  part  of  the  analog  computer  capability  for  other  uses,  it  was  decided  to 
differentiate  the  functions  of  the  two  experimental  systems  in  the  above  manner 
in  general  to  provide  a more  readily  comprehensible  approach  to  a complex  prob- 
lem and  specifically  to  separate  out  the  simulated  AV-8A  control  input  from  the 
actual  X-22A  control  input  for  the  purpose  of  control  usage  measurements.  This 
separation  of  functions  between  the  VSS  and  analog  computer  is  illustrated 
schematically  in  Figure  5-1. 

The  aerodynamic/propulsive  force  and  moment  characteristics  of  che 
AV-8A  in  transition  were  mechanized  on  the  VSS  in  terms  of  perturbations  from 
a reference  transition  rather  than  by  applying  the  fixed  operating  point  simula- 
tion technique  which  involves  perturbations  of  the  motion  and  control  variables 
from  initial  trim  conditions.  The  perturbation  technique  involves  not  only 
feedback  gains  and  control  gearings  which  vary  as  a function  of  flight  con- 
dition (Section  IV)  but  also  motion  feedback  and  control  variables  expressed 
in  terms  of  perturbations  from  reference  values  which  are  also  functions  of 
flight  condition.  As  applied  to  the  AV-8A  simulation  problem,  the  technique 
ideally  yields  a valid  simulation  of  the  aircraft  dynamics  for  "small"  pertur- 
bations of  the  motion  and  control  variables  about  their  reference  values  at 
various  fixed  operating  points  along  the  reference  transition  as  well  as  a 
faithful  reproduction  of  the  AV-8A  trim  changes  with  flight  condition.  For 
example,  with  the  VSS  engaged,  the  actual  X-22A  pitch  control  input  ( ) 

may  be  expressed  as  (see  Figure  5-1): 
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Figure  5-1  AV-8A  MODEL  MECHANIZATION 
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where  Ae  (A)  = actual  X-22A  "trim"  pitch  control  position  along  the  refer- 
ence transition;  a function  of  A (Figure  5-2). 
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*4/^  : reference  values  of  u , -as  as  func- 
tions of  A . (Figure  5-3). 


perturbations  from  AV-8A  "trim"  control 
positions  along  reference  transition 
(Figure  5-4) . 


VSS  gains;  functions  of  A (Table  4-5b). 


A similar  expression  may  also  be  written  for  the  actual  X-22A  thrust  control 
input  (.Acs  )•  A solution  of  the  pitch  and  thrust  control  equations  with  the 
X-22A  in  "trim"  (i.e.,  Acs  = A e ; Acg  - Ac  ) along  the  reference  transition 
(Au.  = Aur  - (J.  = 0)  yields  = A6r  = 0 with  the  restriction  that: 


for  any  value  of  A 


Hence  for  the  values  of  longitudinal  control  gearings  mechanized  (Table  4-8), 
the  "trim"  positions  of  the  evaluation  pilot's  control  without  SCAS  inputs  are 
always  the  simulated  AV-8A  trim  positions  programmed  on  the  analog  computer; 
the  presence  of  SCAS  inputs  will  alter  the  trim  control  positions  in  a manner 
similar  to  that  which  would  occur  on  the  actual  aircraft. 


As  discussed  in  general  in  Section  IV  and  specifically  above,  the 
gain  schedulings  required  for  the  AV-8A  simulation  are  performed  as  functions 
of  X-22A  duct  angle,  hence  airspeed.  However,  the  AV-8A  throttle  control  deriv- 
atives ( r ) are  in  general  also  functions  of  AV-8A  nozzle  angle 

(8j  );  consequently,  for  a more  accurate  simulation  of  the  aircraft  response  to 
throttle  inputs,  an  additional  function  generator  may  be  required  to  alter  the 
magnitude  of  the  AS  r signal  (Figure  5-1)  to  the  VSS  and  collective  elevon  logic 
as  a function  of  6j  in  order  to  achieve  the  required  variation  in  throttle  deriv- 
atives with  AV-8A  configuration  change. 
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Section  IV  also  discusses  the  requirement  for  a third  controller  in 
the  longitudinal  degrees  of  freedom.  This  controller,  the  collective  use  of 
the  elevons,  has  been  incorporated  in  the  simulation  facility  by  an  addition  to 
the  simulator  digital  computer  program  to  account  for  the  effects  of  collective 
elevon  deflection  on  the  basic  X-22A  and  the  programming  of  the  logic  to  drive 
the  elevons  in  a collective  fashion  on  the  analog  computer.  The  alteration  of 
the  digital  computer  program  consists  simply  of  the  addition  of  a single  term 
each  to  the  expressions  for  the  X and  £ aerodynamic/propulsive  forces: 

and  ^ec  respectively;  the  pitching  moment  due  to  col- 

lective elevon  deflections  is  assumed  to  be  negligible.  The  values  of  the 
collective  elevon  effectiveness  derivatives,  estimates  based  upon  X-22A  wind 
tunnel  data  and  geometry,  are  programmed  in  tabular  form  similar  to  that  used 
for  the  other  stability/control  parameters;  selected  values  of  these  derivatives 
are  presented  in  Table  5-2. 


TABLE  5-2 

COLLECTIVE  ELEVON  EFFECTIVENESS  DERIVATIVES 


A 

Deg 

X /> 

&£L 

Ft/Sec2/Deg 

Ft/Sec^/Deg 

0 

0.0 

-0.4 

15 

-0.1 

-0.39 

30 

-0.2 

-0.35 

45 

-0.28 

-0.28 

65 

-0.36 

1 

o 

90 

-0.40 

0.0 

NOTE:  [Positive  Down  ( A = 0) 

ec  [Positive  Forward  ( A = 90) 

The  logic  to  drive  the  elevons  collectively  is  currently  programmed  on  the 
analog  computer  as : 

*«.  - a.sT 

where  an<*  are  functions  of  A as  presented  in  Section  4.2.4 

and  AeL  is  the  collective  elevon  input  signal  to  the  digital  computer; 
currently  limited  to  +4.5  deg. 
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5.3 


SIMULATOR  CHECKOUT 


Both  static  and  dynamic  tests  were  conducted  to  ensure  the  proper 
operation  of  the  various  elements  of  the  simulator  facility  when  integrated 
as  indicated  in  Figure  5-1  to  provide  a simulation  of  the  AV-8A  aircraft.  The 
static  tests  consisted  primarily  of  operational  checks  of  the  analog  computer 
and  VSS  function  generators  by  inserting  known  values  of  the  state  and  control 
variables,  measuring  the  appropriate  output,  and  comparing  this  value  to  the 
desired  value  at  each  duct  angle;  the  operation  of  the  collective  elevons  was 
checked  out  using  a similar  input-output  technique.  Dynamic  testing  consisted 
of  both  fixed  operating  point  time  history  comparisons  and  piloted  decelerat- 
ing approaches.  Longitudinal  and  lateral-directional  responses  to  control 
inputs  at  trimmed  velocities  of  0,  60,  and  120  knots  were  recorded  and  com- 
pared to  the  time  histories  presented  in  Figures  4-1  and  4-4;  the  dynamic 
response  characteristics  of  the  simulator  were  judged  to  match  the  desired 
characteristics  satisfactorily.  The  final  test  of  the  AV-8A  simulation  con- 
sisted of  simulated  decelerating  approaches  with  the  basic  AV-8A  plus  SAS  (Ap- 
pendix II)  along  the  Task  III  approach  profile  using  the  automatic  duct  rota- 
tion feature  of  the  X-22A,  the  control  director  electronic  display  format 
(ED-3)  with  the  RATE  system  director  logic,  and  the  throttle  control  of  thrust 
magnitude  (Reference  31);  no  obvious  system-related  problems  were  revealed. 


« 
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Section  6 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  general  conclusion  from  this  study  is  that  it  is  feasible  to 
simulate  in  the  X-22A  aircraft  an  AV-8A  class  VTOL  performing  terminal  area 
operations.  As  has  been  described  in  this  report,  most  of  the  methodology 
for  performing  such  a simulation  has  been  developed  and  mechanized  for  fur- 
ther experimental  investigation  on  the  X-22A  ground  simulator;  it  has  been 
demonstrated  that  excellent  simulation  fidelity  below  60  kt.  is  possible  and 
that  the  mechanization  of  the  simulation  can  be  performed  with  the  X-22A  equip- 
ment . 


Specifically,  the  study  objectives  were  achieved  by  the  accomplish- 
ment of  the  following  tasks: 

• A generic  model  of  the  stability  and  control  derivatives  of  the 
AV-8A  was  developed,  primarily  from  the  wind  tunnel  data  given 
in  Reference  2.  The  model  was  developed  for  an  angle  of  attack 
of  8 degrees  and  glide  slope  angle  of  5 degrees  to  be  consistent 
with  previous  simulations  (References  5,  6);  values  of  the  deri- 
vatives were  computed  for  six  trimmed  velocities  ( l/a  =0,  30, 

60,  80,  100,  120  kt).  The  model  is  linear  around  a set  of  de- 
veloped reference  conditions  which  are  a function  of  trim  velocity 
only.  Comparisons  of  the  characteristic  roots  with  available 
data  in  Reference  3 indicate  excellent  agreement  at  all  six  trim 
velocities . 

• The  variable  stability  gains  for  the  lateral-directional  simu- 
lation were  computed  using  a Calspan-developed  procedure  which 
matches  desired  modal  characteristics  (Reference  28).  Both  feed- 
back and  feedforward  (gearing)  gains  to  the  lateral  stick  and 
rudder  pedals  were  calculated  for  all  six  flight  conditions. 
Excellent  matching  of  time  history  responses  in  sideslip,  roll 
rate,  roll  attitude,  and  yaw  rate  to  both  aileron  and  rudder 
inputs  was  achieved  at  all  flight  conditions.  The  only  discrep- 
ancy is  the  lateral  acceleration  response,  which  is  signifi- 
cantly higher  in  the  simulation  than  in  the  AV-8A  model;  this 
difference  arises  because  of  the  larger  Ya  in  the  X-22A  which 
cannot  be  altered  with  the  variable  stability  system.  Flying 
qualities  parameters  of  the  AV-8A  model  and  the  simulation  were 
computed  for  comparison  with  M1L-F-83300  (Reference  23);  for  all 
flight  conditions,  these  parameters  were  either  exactly  the 

same  or  indicated  equivalent  levels  of  flying  qualities. 

• Several  methods  of  computing  the  variable  stability  gains  for 
the  longitudinal  simulation  were  examined.  An  extension  of 
model-in-the-performance-index  optimal  control  theory  to 
allow  computations  with  unstable  models  was  developed  and 
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selected  as  the  most  useful  method  for  computing  the  feedback 
gains  assuming  two  controllers  (longitudinal  stick  and  col- 
lective blade  pitch).  The  major  criterion  in  the  iterative 
application  of  this  procedure  to  compute  feedback  gains  was 
that  the  characteristic  roots  of  the  simulation  match  those  of 
the  model;  this  criterion  was  met  well  at  all  six  velocities. 

The  gearing  (feedforward)  gains  were  computed  using  an  equa- 
tion-error least-squares  technique  that  reproduces  the  model 
control  matrix  (and  hence  the  initial  accelerations  to  control 
inputs)  exactly  if  three  controllers  are  used.  For  these  com- 
putations, it  was  necessary  to  assume  the  expansion  of  the 
variable  stability  system  capability  to  include  collective  ele- 
vons,  thereby  providing  a third  controller.  This  expansion  is 
necessary  because  of  the  markedly  different  thrust  inclinations 
at  a given  speed  between  the  X-22A  and  AV-8A  model:  to  match 

initial  responses  to  simulated  thrust  inputs  at  the  higher 
speeds  (^  = 100,  120),  which  are  low  duct  angles  for  the  X-22A, 
the  collective  elevons  are  required  to  match  the  simulated 
thrust  inclination.  As  a result  of  using  three  controllers  for 
the  gearing  calculations,  the  simulated  control  matrix  is  an 
exact  duplicate  of  the  AV-8A  model  control  matrix.  Excellent 
matching  of  time  history  responses  in  longitudinal  velocity, 
vertical  velocity,  pitch  rate,  and  pitch  attitude  to  both  simu- 
lated pitch  stick  and  throttle  inputs  was  achieved  for  the  three 
low-speed  flight  conditions  ( =0,  30,  60  kt) . Good  match- 

ing of  the  responses  to  pitch  stick  inputs  is  also  obtained  for 
the  three  higher  speed  conditions  ( l/0  = 80,  100,  120  kt) , but 

vertical  velocity  responses  to  simulated  throttle  inputs  do  not 
match.  The  discrepancy  is  attributable  to  significant  differ- 
ences between  the  simulated  values  and  AV-8A  model  values  of  the 
stability  derivatives  jgu  and  , which  cannot  be  modified  for 
high  speeds  (low  duct  angles)  with  the  two  controllers  used  for 
the  feedbacks.  Although  the  matching  of  the  characteristic  roots 
yields  equivalent  levels  of  flying  qualities  in  terms  of  MIL-F- 
83300  requirements  for  the  simulation  and  the  AV-8A  model  at  the 
higher  speed  conditions,  it  is  recommended  that  additional  stu- 
dies be  performed  to  improve  the  fidelity  of  the  throttle  re- 
sponses for  these  conditions;  the  excellent  matches  to  both  stick 
and  throttle  inputs  for  the  low  speed  cases  imply  good  simulation. 

• The  X-22A  ground  simulator  facility  was  updated  and  modified  to 
reflect  as  closely  as  possible  the  aircraft  characteristics  and 
mechanization.  The  mathematical  model  of  the  basic  X-22A,  as 
programmed  in  the  simulator's  digital  computer,  was  updated  by 
incorporating  values  of  stability  and  control  derivatives  that 
were  identified  from  flight  data  during  the  Task  ill  X-22A  exper- 
iment (Reference  31).  Incorporation  of  the  airborne  analog  com- 
puter into  the  operation  of  the  ground  simulator  was  accom- 
plished; this  computer  performs  the  guidance  and  control  mechan- 
ization functions  in  flight,  and  its  incorporation  in  the 
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simulator  permits  exact  duplication  of  these  functions.  Finally, 
the  ground  simulator's  variable  stability  system  - an  exact 
duplicate  of  that  in  the  aircraft  — was  made  operational,  there- 
by allowing  duplication  of  the  function  generation  for  the  gains 
and  reference  conditions  that  will  be  required  in  flight. 

• The  simulation  of  the  AV-8A  model  was  mechanized  on  the  ground 
simulator.  The  underlying  philosophy  of  the  mechanization  is 
to  perform  the  simulation  of  the  basic  AV-8A  with  the  VSS  as 
much  as  possible,  and  to  use  the  airborne  analog  computer  pri- 
marily to  implement  control  system  designs;  although  some  of 
the  control  system  functions  could  be  performed  on  the  VSS,  the 
separation  was  selected  to  permit  direct  measurement  of  AV-8A 
control  usage  for  the  various  control  systems  and  to  provide  a 
rational  method  for  dealing  with  a complex  problem.  According- 
ly, the  gain  scheduling  and  state  variable  reference  condition 
variations  were  programmed  on  function  generator  cards  in  the 
VSS  as  a function  of  duct  angle  (and  hence  trim  velocity),  and 
the  operation  of  the  VSS  was  modified  to  (1)  permit  reference 
values  for  longitudinal  and  vertical  velocity  to  be  introduced, 
and  (2)  permit  longitudinal  velocity,  pitch  rate,  and  pitch 
attitude  feedbacks  to  the  collective  blade  control.  Reference 
or  trim  values  of  the  simulated  throttle  and  longitudinal  stick 
are  introduced  at  the  interface  between  the  analog  computer  and 
the  VSS.  The  feel  system  is  used  to  reproduce  the  force- 
displacement  characteristics  of  the  AV-8A  cockpit  controllers;  a 
fore-and-aft  controller  simulates  the  throttle,  and  a smaller 
fore/aft  controller  will  simulate  the  nozzle  angle  control.  For 
checkout  purposes,  a simplified  linear  approximation  to  the 
AV-8A  stability  augmentation  system  is  programmed  on  the  analog 
computer.  Checkouts  of  the  simulation  mechanization  have  been 
conducted  with  satisfactory  results. 

It  is  recommended  that  additional  attention  be  given  during  the  ini- 
tial stages  of  the  Task  IV  ground  simulation  experiment  to  the  following  prob- 
lems discovered  during  this  study: 

• Responses  to  throttle  inputs  at  the  higher  speeds.  It  may  be 
necessary  to  (1)  consider  limited  feedback  to  the  collective 
elevons  or  (2)  trade  off  initial  acceleration  responses  against 
more  closely  matched  velocity  responses  at  a specified  time. 

• Trade  off  sideslip  response  matching,  which  is  excellent  with 
the  current  simulation,  against  lateral  acceleration  matching 
at  the  pilot's  station. 

• Initial  deceleration  for  large  nozzle  angle  changes.  The  cur- 
rent simulation  can  match  the  majority  of  the  AV-8A  decelera- 
tion profile  for  nozzle  changes,  which  is  essentially  exponen- 
tial, by  duct  angle  changes.  An  exception  is  the  initial 
deceleration,  which  may  require  washed-out  collective  blade 
pitch  inputs  as  a function  of  nozzle  change. 
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• Data  deficiencies.  Additional  AV-8A/B  data,  particularly  for 
reference  conditions  and  throttle  effectivenesses,  are  desir- 
able . 
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GLOSSARY  OF  SYMBOLS 


**•  C ) 

J 

J- 

J 

K.  < 

L 

L\  ) 


error  term 

characteristic  matrix,  open  loop 
characteristic  matrix,  closed  loop 
acceleration  due  to  gravity,  32.2  ft/sec^ 
control  matrix 
measurement  matrix 
identity  matrix 

moment  of  inertia  about  body  ( )-axis  (ft-lb/sec“) 

product  of  inertia  in  body  axes  (ft-lb/sec^) 

complex  variable,  l4T 

gearing  matrix 

integrand  (section  4.2.2) 

feedback  matrix 

aerodynamic  moment  about  body  X-axis  (ft-lb) 
dimensional  rolling  moment  derivative 


± ( 

~'\3L  + 

I** 

3 N " 

1*  \ 

rxij 

[3(  ) 

d(  ) . 

A' 

i 


aerodynamic  moment  about  body  Y'-axis  (ft-lb) 

dimensional  pitching  moment  derivative 
1 dM  rad/sec2 

* fy  ~T) 

integer 

body  axis  accelerations  (ft/sec^) 
aerodynamic  moment  about  £ -axis  (ft-lb) 
dimensional  yawing  moment  derivative 


/ - r**  \ 

[*( ) 

j[K  3(  ) 

numerator  of  the  l/j  (s)  transfer  function 
body-axis  roll  rate  (deg/sec,  rad/sec) 
steady-state  augmented  control  gain  matrix 
body-axis  pitch  rate  (deg/sec,  rad/sec) 
state  weighting  matrix 
body-axis  yaw  rate  (deg/sec,  rad/sec) 
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GLOSSARY  OR  SYMBOLS  (Cont'd) 


AC  ) 
A Cs  ) 


control  weighting  matrix 
Laplace  operator  cr  £ j cO 
time  (sec) 

velocity  along  body  X-axis  (ft/sec) 
control  vector 

velocity  along  body  V-axis  (ft/sec) 

velocity  (kt,  ft/sec) 

velocity  along  body  2 -axis  (ft/sec) 

aerodynamic  force  along  bodyX.Y',  2 axis,  respectively  (11 
2(  ^dimensional  longitudinal,  lateral,  vertical  force  derivati\ 
1 t?X,  Yt  or  2 f ft/sec2  ^ 
m 9 ( ) \ (“)  / 

state  vector 

angle  of  attack  (deg,  rad) 
angle  of  sideslip  (deg,  rad) 
flight  path  angle  (deg) 
pilot's  controller  position 

e,£,es,£S  longitudinal  stick  (in.),  positive  aft 

a,.  A, as.  AS  lateral  stick  (in.),  positive  right 

r,  rp,  RP  rudder  pedals  (in.),  positive  right 

T,  cs,  CS  throttle  or  collective  stick  (in.),  positive 
increased  thrust 

perturbation  term  ( ) - ( )0,  units  of  ( ) 
characteristic  equation 

X-22A  safety  pilot's  or  VSS  controller  position 

e,£,Cs,  ES  longitudinal  (in.),  positive  nose  up 

<Z.A,as,A5  lateral  (in.),  positive  right  wing  down 

r , rP>  RP  directional  (in.),  positive  nose  right 

CS,  CS  thrust  (in.),  positive  increased  thrust 

cl.£L  collective  elevons  (deg),  positive  trailing  e 

down  and/or  forward 

damping  ratio 

damping  ratio  of  dutch  roll  characteristic  roots 
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Jl 

cr 
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u> 
^cL 
00 n 

“V 
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v« 


damping  ratio  of  numerator  roots  in  <P/8as  transfer  function 

pitch  attitude  (deg,  rad) 

nozzle  angle  Mjri  engine  axis  (deg) 

nozzle  angle  Mrrt  fuselage  reference  line  (deg) 

duct  angle,  measured  from  horizontal  (deg) 

adjoint  variable 

eigenvalue  matrix  of  ( ) 

real  part  of  Laplace  operator 

generalized  time  constant 

roll  mode  time  constant  (sec) 

spiral  mode  time  constant  (sec) 

Euler  roll  attitude  (deg,  rad) 
magnitude  of  ratio  for  Dutch  roll  mode 

imaginary  part  of  Laplace  operator  (rad/sec) 
undamped  natural  frequency  of  Dutch  roll  mode  (rad/sec) 
undamped  natural  frequency  (rad/sec) 

undamped  natural  frequency  of  numerator  roots  in  / 8ag 

transfer  function  (rad/sec) 

partial  derivative  of  L with  respect  to  j 

gradient  (vector)  with  respect  to  U. 


Subscripts 


( 

>xs 

least-squares  solution 

( 

model 

( 

optimal  control  solution 

( 

) R 

reference  trajectory  value 

( 

^Lfv 

lateral -directional 

( 

^ L OA/S 

longitudinal 

( 

V 

plant 

( 

reference  value  (also  trim) 

( 

^ trim 

trim  value 

( 

desired  steady-state  responses  (section  4.2.4) 
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Superscripts 

C) 

( )T 

V-/ 


( ) 

A 

( 


( > 
( ) 


derivative  with  respect  to  time 
matrix  transpose 
matrix  inverse 
combined  matrices 

transformed  problem  (section  4.2.3) 
value  of  ( ) is  degrees 
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Appendix  I 

BASIC  AV-8A  MODEL  STABILITY/CONTROL  DERIVATIVES 


The  following  plots  of  both  longitudinal  and  lateral-directional 
stability  and  control  derivatives  as  functions  of  airspeed  represent  the  re- 
sults of  the  AV-8A  data  gathering  and  analysis  portions  of  the  preliminary 
program.  The  AV-8A  mathematical  model  presented  in  Section  2 is  based  pri- 
marily upon  these  data.  An  explanation  of  the  symbology  used  to  denote  the 
source  of  each  data  point  follows: 

O NRC,  NAE  LR-500:  National  Research  Council  of  Canada  simula- 

tion of  PI 127  (V  = 50,  80  kts) 

□ MDC  A1410:  MCAIR  AV-8A  aerodynamic  data,  slow  approach 

( Oj  = 70°) 

()  MDC  A1410:  MCAIR  AV-8A  aerodynamic  data,  decelerating  tran- 
sition = 90  ) 

A NR  72H-268:  NR  simulation  of  AV-8A  "validated"  at  V = 25  knots 

(indicated  by  shaded  symbol) 

^ MDC  A3618:  MCAIR  Landing  Aids  study;  revised  estimates  of 

O MDC  A3580:  MCAIR  AV-8A  model  for  AFCS  design;  estimates  of 

A«/’  ' Lv  ’ /Vv'  Presented  here. 
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Appendix  II 

AV-8A  SAS  CHARACTERISTICS 


General 

The  AV-8A  SAS  is  a limited  authority  three-axis  system  consisting  basically 
of  pitch  and  roll  rate  damping  augmentation  and  a turn  coordination  feature  in 
the  yaw  axis.  The  pitch  SAS  operates  through  the  stabilator  and  the  downward- 
blowing tail  reaction  control  system  (RCS)  valve  only;  this  characteristic 
introduces  an  anomaly  in  the  longitudinal  SAS-on  dynamics  in  hover  and  low  speed 
flight  which  will  be  discussed  later  in  this  appendix.  Both  the  ailerons  and  the 
roll  RCS  valves  are  operated  by  the  roll  SAS;  however,  the  yaw  SAS  uses  only 
the  dual  yaw  RCS  valve  at  the  tail,  thus  isolating  the  pilot's  rudder  pedals 
from  SAS  inputs.  Apparently  the  pitch  and  roll  SAS  inputs  also  move  the  pilot's 
control  stick.  The  details  of  the  SAS  dynamic  characteristics  are  discussed 
in  the  following  sections. 

Pitch  SAS 


The  AV-8A  pitch  SAS  consists  of  the  feedback  of  body-axis  pitch  rate  ( o ) 
to  the  stabilator  and  aft  pitch  RCS  valve.  This  feedback  signal  is  limited  to 
a value  which  corresponds  to  ±1.5  deg  of  stabilator  travel  ( & H ) or  approxi- 
mately *0.84  in.  of  longitudinal  stick  travel  ( Ses  ).  The  stabilator  angle- 
to-pitch-rate  transfer  function  is  as  follows: 
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During  a decelerating  transition  or  slow  approach  the  dominant  pitch 
control  moment  is  supplied  by  the  RCS.  Since  the  pitch  SAS  operates  through  the 
stabilator  which  is  linked  directly  to  the  aft  pitch  RCS  valve,  the  stabilizing 
effect  of  the  SAS  inputs  is  significantly  reduced  when  the  stabilator  setting 
is  such  that  the  valve  is  closed  (Reference  21).  Unfortunately  this  phenomenon 
occurs  during  the  hover  and  low  speed  portion  of  the  landing  approach.  Figure  A-l 
demonstrates  that  for  airspeeds  below  55  knots  the  trim  value  of  , hence 

6fs  * is  such  that  the  rear  RCS  valve  is  closed  and  the  front  RCS  valve 
open.  As  the  transition  progresses,  larger  and  larger  error  signals  are  re- 
quired to  open  the  rear  valve;  in  fact,  for  airspeeds  less  than  ~5  knots 
the  maximum  value  of  the  limited  error  signal  (+0.84  in.)  will  not  open  the 
valve  at  all  without  a pilot  control  input. 

Therefore  for  airspeeds  less  than  ~55  knots  the  pitch  SAS  provides  little 
or  no  restoring  moment  in  response  to  a nose-down  pitch  rate  (stabilator  effect 
only)  and  provides  a significant  restoring  moment  (stabilator  plus  RCS)  only 
in  response  to  nose-up  pitch  rates  greater  than  some  threshold  value  which 
increases  with  decreasing  airspeed. 


Roll  SAS 

The  roll  SAS  consists  of  the  feedback  of  body-axis  roll  rate  ( f ) to  the 
ailerons  and  roll  RCS  valves.  This  feedback  signal  is  limited  to  a value  which 
corresponds  to  ± 2 deg  of  aileron  travel  ( SA  ) or  approximately  +0.68  in.  of 
lateral  stick  travel  ( Sa_s  ).  Hie  aileron-to-roll-rate  transfer  function  is 
as  follows: 


-1/93 
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Since  ^as/^A  * 0*34  in. /deg, 
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The  roll  RCS  valves  are  linked  to  the  ailerons  so  that  at  deg  of  aileron 
deflection  the  downward-blowing  valve  is  fully  open.  During  the  remainder  of 
the  total  aileron  travel  (±  12.0  deg)  the  upward -blowing  valve  on  the  opposite 
wing  opens.  The  upward-blowing  valve  is  less  effective  in  providing  a rolling 
moment;  hence  the  functional  relationship  of  rolling  moment  (L)  and  aileron 
deflection  is  "kinked"  at  SA  = +7°;  i.e.. 
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Yaw  SAS 

The  yaw  SAS  consists  of  the  feedback  of  lateral  acceleration  as  measured 
at  the  accelerometer  ( ) and  washed-out  yaw  rate  ( ) , to  the  yaw 

RCS  valve  and  an  aileron-to-yaw  RCS  valve  interconnect.  The  sum  of  the  feed- 
back and  feedforward  signals  is  limited  to  a value  which  corresponds  to  ±50% 
yaw  RCS  valve  opening  (±5  deg  equivalent  rudder  ( 6*f  ) ) or  approximately 

*T0.7  in  rudder  pedal  travel  ( igp  ).  The  equivalent  rudder-to-lateral- 
acceleration  transfer  function  is  as  follows: 
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The  equivalent -rudder-to-yaw-rate  transfer  function  is: 


or 
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Finally,  the  equivalent-rudder-to-aileron  transfer  function  is: 


or 


-0.58 


(T%r)(^b)  • 

F/L  Ten  HCrUATOR. 


deg 

deg 


- 0.24.  ( 


)r 


)■ 


in. 

in. 


95 


Pitch  RCS  Schedule  • Transition  Trim  Stick  Position 
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